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An abundance of astrophysical evidence indicates that the bulk of matter in the universe is made up of massive, electrically 
neutral particles that form the dark matter (DM) . While the density of DM has been precisely measured, the identity of the 
DM particle (or particles) is a complete mystery. In fact, within the laws of physics as we know them (the Standard Model, 
or SM), none of the particles have the right properties to make up DM. Remarkably, many new physics extensions of the 
SM - designed to address theoretical issues with the electroweak symmetry breaking sector - require the introduction of new 
particles, some of which are excellent DM candidates. As the LHC era begins, there are high hopes that DM particles, along 
with their associated new matter states, will be produced in pp collisions. We discuss how LHC experiments, along with other 
DM searches, may serve to determine the identity of DM particles and elucidate the associated physics. Most of our discussion 
centers around theories with weak-scale supersymmetry, and allows for several different DM candidate particles. 



1. Introduction 

The LHC program has been described as the great- 
est experiment ever to be mounted in physics. Cer- 
tainly this seems to be true on many different levels: 
the largest, costliest, most massive detectors; the most 
collaborators per experiment; the highest energy reach 
of any accelerator experiment. The intellectual stakes 
of the LHC program are enormous: on the theory side, 
the extreme sensitivity of the scalar sector of the Stan- 
dard Model (SM) to very high scale physics beckons 
for new physics at the weak scale (~ 100 — 1000 GeV), 
possibly ushering in a new paradigm for the laws of 
physics. 

We discuss how LHC experiments may serve to val- 
idate the extended Copernican principle. In previous 
times, we have learned that the earth is not the center 
of the solar system, that our galaxy is not the entire 
universe, and that we do not live in any special place 
or time. Now, due to an impressive accumulation of as- 
trophysical data, we learn that our star, our planet, and 
ourselves are not even made up of the dominant form 
of matter in the universe. It now appears that most of 
the matter in the universe - the so-called dark matter 
(DM) - must consist of massive, electrically and (likely) 
color neutral particles that were produced with non- 
relativistic velocities (cold DM or CDM) in the early 
universe. None of the particles of the SM have the 
right properties to make up CDM. Thus, CDM consti- 
tutes decisive evidence for physics beyond the Standard 
Model [11]! 

Compelling arguments suggest the CDM particle is 
linked to the weak nuclear interactions, and further. 



that it has a mass of order the weak scale: ^ 100—1000 
GeV. This is often referred to as the WIMP miracle, 
and the dark matter particles referred to as WIMPS 
(weakly interacting mass particles). Many attractive 
theoretical scenarios designed to ameliorate the ex- 
treme sensitivity of the scalar sector of the SM to radia- 
tive corrections, naturally include candidates for CDM 
particles with weak scale masses that interact with or- 
dinary matter with cross sections comparable to those 
for weak nuclear interactions. Regardless of its origin, 
if CDM is composed of WIMPs, then it may be pos- 
sible to produce and study the DM particle (s) directly 
at the LHC. In fact, the LHC may well turn out to be 
a DM factory, where the nature of DM particles and 
their properties might be studied in a controlled envi- 
ronment. In any collider experiment, WIMPS would 
be like neutrinos in that they would escape the detec- 
tor without depositing any energy in the experimental 
apparatus, resulting in an apparent imbalance of energy 
and momentum in collider events. While WIMPs would 
manifest themselves only as missing (transverse) energy 
at collider experiments, it should nevertheless be pos- 
sible to study the visible particles produced in WIMP- 
related production and decay processes to study the 
new physics associated with the WIMP-sector. 

Indeed, there exists a real possibility that much of 
the mystery surrounding DM and its properties can be 
cleared up in the next decade by a variety of exper- 
iments already operating or soon-to-be deployed. In 
this effort, experiments at the LHC will play a crucial 
role. There are - in tandem with LHC- a variety of 
other dark matter search experiments already in oper- 
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ation, or in a deployment or planning phase. Direct 
Detection (DD) experiments seek to directly measure 
relic DM particles left over from early stages of the 
Big Bang. These DD experiments range from terres- 
trial microwave cavities that search for axions via their 
conversion to photons, to crystaline or noble liquid tar- 
gets located deep underground that search for WIMP- 
nucleon collisions. 

DM can also be searched for in indirect detection 
(ID) experiments. In ID experiments, one searches 
for WIMP- WIMP annihilation into various SM parti- 
cles including neutrinos, gamma rays and anti-matter. 
Clearly, this technique applies only if the DM is self- 
conjugate, or if DM particles and anti-particles are 
roughly equally abundant. One ID search method in- 
volves the use of neutrino telescopes mounted deep un- 
der water or in polar ice. The idea is that if relic 
WIMPs are the DM in our galactic halo, the sun (or 
earth) will sweep them up as they traverse their galactic 
orbits, and gravitationally trap these in the central core 
where they can accumulate, essentially at rest, to den- 
sities much higher than in the Milky Way halo. These 
accumulated WIMPS can then annihilate one with an- 
other into SM particles with energies E ~ towimp- 
Most SM particles would be immediately absorbed by 
the solar material. However, neutrinos can easily es- 
cape the sun. Thus, WIMP annihilation in the sun re- 
sults in an isotropic flux of high energy neutrinos from 
the solar core ~ these energies are impossible to produce 
via conventional nuclear reactions in the sun - some of 
which would make it to earth. These neutrinos ocas- 
sionally interact with nuclei in ocean water or ice and 
convert to a high energy muon, which could then be de- 
tected via Cerenkov radiation by photomultiplier tubes 
that are parts of neutrino telescopes located within the 
medium. 

Another possibility for ID is to search for the by- 
products of WIMP annihilation in various regions of 
our galactic halo. Even though the halo number den- 
sity of WIMPs would be quite low, the volume of the 
galaxy is large. Ocassionally one expects relic WIMP- 
WIMP annihilation to SM particles. The trick is then 
to look for rare anti-matter production or high energy 
gamma ray production from these WIMP halo annihila- 
tions. A variety of land-based, high altitude and space- 
based anti-matter and gamma ray detectors have been 
or are being deployed. The space-based Pamela exper- 
iment is searching for positrons and anti-protons. The 
land-based HESS telescope will soon be joined by the 
GLAST satellite in the search for high energy gamma 
rays. While high energy anti-particles would provide a 
striking signal, these lose energy upon deflection when 



traversing the complicated galactic magnetic field, and 
so can only be detected over limited distances. Gamma 
rays, on the other hand, are undeflected by magnetic 
fields, and so have an enormous range. Moreover, these 
would point back to their point of origin. Thus, the 
galactic center, where dark matter is expected to ac- 
cumulate at a high density, might be a good source of 
GeV-scale gamma rays resulting from WIMP-WIMP 
annihilation to vector boson {V — W, Z) pairs or to 
quark jets, followed by {V —>-)q — > tt*' ^ 77 after 
hadronization and decay. 

If WIMPs and their associated particles are discov- 
ered at the LHC and/or at DD or ID search experi- 
ments, it will be a revolutionary discovery. But it will 
only be the beginning of the story as it will usher in 
a new era of dark matter astronomyl The next logical 
step would be the construction of an e+e~ collider of 
sufficient energy so that WIMP (and related particles) 
can be produced and studied with high precision in a 
clean, well-controlled experimental environment. The 
precise determination of particle physics quantities as- 
sociated with WIMP physics will allow us to deduce the 
relic density of these WIMPS within the standard Big 
Bang cosmology. If this turns out to be in agreement 
with the measured relic density, we would have direct 
evidence that DM consists of a single component. If 
the predicted relic density is too small, it could make 
the case for multiple components in the DM sector. If 
the predicted density is too large, we would be forced 
to abandon the simplest picture and seek more com- 
plicated (non-thermal) mechanisms to account for the 
measurement, or deduce that this detected WIMP itself 
is unstable. The determination of the properties of the 
DM sector will also serve as a tool for a detailed mea- 
surement of astrophysical quantities such as the galac- 
tic and local WIMP density and local velocity profiles, 
which could shed light on the formation of galaxies and 
on the evolution of the universe. 

2. Evidence for dark matter 

Dark matter in the universe was first proposed in the 
1930s by astronomer Fritz Zwicky['2]. In the 1970s and 
on, evidence for DM accrued at an accelerating pace. 
Here we discuss the major classes of evidence for DM 
in the universe. 

• Galactic clusters: In the 1930s, Zwicky studied 
nearby clusters of galaxies, bound to each other 
by gravity in spite of the expansion of the uni- 
verse. Using arguments based on the virial the- 
orem from classical mechanics, Zwicky concluded 
there was not enough visible mass within the 
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galactic clusters to successfully bind them; he 
thus concluded that there must be large amounts 
of non-luminous, or dark matter, existing within 
the clusters. 

• Rotation curves: In the 1970s, V.C. Rubin and 
W.K. Ford[ 13] began an intensive study of the 
rotation curves of galaxies. They were able to 
measure stellar velocity as a function of distance 
from the galactic center. With most of the visi- 
ble matter concentrated in or around the galactic 
center, one expects the stellar rotational velocities 
to fall off with distance from the galactic center 
in accord with Newtonian gravitation. Instead, 
the stellar velocities tended to flatness out to the 
furthest distances which could be probed. This is 
in accord with a diffuse halo of dark particles sur- 
rounding the galaxy out to the furthest distances. 

• Lensing: In General Relativity, the path of light 
through space-time is bent, or "lensed" as it 
passes by a large mass distribution. Lensing ef- 
fects are observed when light from distant galax- 
ies or clusters passes by large mass distributions. 
Numerous studies of both strong and weak (sta- 
tistical) lensing show the presence of large quan- 
tities of DM in the universe. 

• Hot gas in clusters: Hot gas bound to clusters of 
galaxies can be mapped out by the emitted x-rays. 
The visible mass in these galaxies would not have 
enough gravity to bind the hot gas, which requires 
additional binding from putative DM. 

• Cosmic microwave background (CMB): Detailed 
studies of anisotropics in the cosmic microwave 
background has resulted a very precisely mea- 
sured CMB power spectrum. The peaks and val- 
leys in this spectrum are extremely sensitive to 
the composition of the universe, and indepen- 
dently show that the universe is comprised of 
about 70% dark energy (DE), 25% DM and 4% 
baryons, along with tiny fractions of neutrinos 
and photons. Thus the "known stuff" makes up 
just about 5% of the content of our Universe. 

• Large scale structure: Measurements of large scale 
structure, when compared to simulations of the 
evolution of structure in the universe, match very 
well with a universe composed of both cold dark 
matter (possibly with some warm DM) and DE. 

• Big Bang nucleosynthesis: One of the triumphs of 
Big Bang cosmology is that given an initially hot. 



dense universe, one can calculate the abundances 
of the light elements produced via nucleosynthesis 
during the first few hundred seconds. The mea- 
sured abundances agree with observation if the 
baryon-to-photon ratio rjB = ^^s/'^7 ~ 6 x 10"^". 
The photon number density is known from ther- 
modynamics, so this implies a baryonic mass den- 
sity of the universe of about ~ 4%, consistent 
with the value independently obtained from CMB 
data discussed above. 

• Distant supernovae probes: Probes of distant su- 
pernovae[ |4j have allowed an extension of the 
Hubble diagram out to redshifts of 2 ^ 1. A 
best fit match to the Hubble diagram indicates 
the presence of both dark energy and dark mat- 
ter in the universe. 

• Colliding galactic clusters: Observation of collid- 
ing clusters of galaxies - a recent example comes 
from the so-called bullet cluster - shows an actual 
separation of dark matter (deduced from lensing) 
from the gaseous halo made of baryonic matter. 
This is exactly what is expected if a vast halo 
of non-interacting dark matter accompanies the 
luminous matter and gas in galactic clusters. 
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Figure 1. Measurements from CMB, large 

scale structure and supernovae plotted in 
the Qa vs. Qmatter plane. Adapted from 

h.ttp://www. astro. washington.edu/astro323/WebLectures/\ 

The KCDM universe: Collating all the data to- 
gether, especially that from CMB, red shifts of high- 
z supernovae, and large scale structure, allows one to 
fit to the composition of the universe. We see from 
Fig. [1] that these very diverse data find consistency 
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amongst themselves, leading to the so-called "concor- 
dance" model for the universe, the KCDM model. 
(Here, A stands for Einstein's cosniological constant, 
which may be the source of the DE). In the KCDM 
model, the universe is composed of about 70% DE, 25% 
DM, 4% baryons with a tiny fraction of neutrinos and 
radiation. The measured abundance of CDM in our 
universe [[5], 

r!ci3M/i' = 0.111^[;;0ii (2a), (1) 

where Vic dm = Pcdm/Pc with pcDM the CDM mass 
density, pc the critical closure density and h is the scaled 
Hubble parameter, serves as a severe constraint on all 
particle physics theories that include a dark matter can- 
didate. Since DM may well consist of more than one 
component, strictly speaking the relic density serves as 
an upper bound Q.xh?' < 0.122 on the density of any 
single component X. We now turn to a discussion of 
some of the particle physics candidates for the DM par- 
ticle X. 
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Figure 2. Dark matter candidates in the mass versus in- 
teraction strength plane, taken from Ref.f^. 



3. DM candidates 

While the evidence for the existence of DM in the 
universe is now very convincing, and while the density 
of dark matter in the universe is becoming precisely 
known, the identity of the dark matter particle(s) is 
a complete mystery. None of the particles in the Stan- 
dard Model have the right properties to make up CDM. 
Many candidates, however, have been proposed in the 
theoretical literature. To appreciate the variety of can- 
didate particles proposed, we list a number of possibil- 
ities. The range of masses and interaction strengths of 
mai»yM^tfefi&ftig&n]dIi^yiW;i&^tJMbSnaFigw[2kkly inter- 
acting neutral massive particles and so are nat- 
ural candidates for the DM in the universe [ [S]. 
It is now known that the usual active neutrinos 
are so light that they could not give rise to the 
observed structure in the Universe because these 
would move faster than the typical galactic escape 
velocity, and so cannot cause the clumping that 
large scale structure simulations require. They 
are usually referred to as hot DM, or HDM, and 
are likely to be a subdominant component of the 
DM in the Universe. There are, however, pro- 
posals for much heavier, cold dark matter gauge 
singlet neutrinos that are not part of the Stan- 
dard Model[[8]. 

• Planck mass black hole remnants: It is possible 
many tiny black holes (BHs) were produced in the 
early universe. Ordinarily, these BHs would de- 
cay via Hawking radiation. However, it has been 



suggested that once they reach the Planck mass, 
quantum gravity effects forbid further radiation, 
making them stable, and hence good CDM can- 
didates [f9]. 

• Q-balls: These objects are topological solitons 
that occur in quantum field theorv[ [TOl lllj . 

• Wimpzillas: These very massive beasts were pro- 
posed to show that viable DM candidates could 
have masses far beyond the weak scale [fT2]. 

• Axions: The symmetries of the QCD Lagrangian 
allow the term - £ 9 ^^^Pt'^F^^- which gives 
rise to CP violation in the strong interactions. 
However, measurements of the neutron electric 
dipole moment (EDM) require Oqcd ^ 10~^°. 
Why this parameter is so much smaller than its 
natural value of ^ 1 is referred to as the strong 
CP problem. The most compelling solution 
to the strong CP problem - the Peccei-Quinn- 
Weinberg-Wilczek solution [ [13] - effectively re- 
places the parameter Oqcd by a quantum field, 
and the potential energy allows the field to re- 
lax to near zero strength. However, a remnant 
of this procedure is that a physical pseudoscalar 
boson - the axion a - remains in the spectrum. 
The axion is an excellent candidate for CDM in 
the universe [ 114] . Its favored mass range is nia ~ 
10~^ — 10"'^ eV, where the lower bound gives too 
high a relic density, and the upper bound comes 



Dark matter and the LHC 



5 



from limits on stellar cooling. Axions have a very 
weak but possibly observable coupling to two pho- 
tons. They are at present being searched for in 
terrestrial microwave cavity experiments such as 
ADMX[fr51. Since they have little direct impact 
on LHC physics, we will not dwell on them in as 
much detail as some other possible candidates. 

WIMPs and the WIMP miracle: Weakly interact- 
ing neutral, massive particles occur in many par- 
ticle physics models where the SM is extended to 
address the physics associated with electroweak 
symmetry breaking (EWSB). If the associated 
new particles sector has a conserved "parity-like" 
quantum number that distinguishes it from the 
SM sector, the lightest particle in this new sec- 
tor is stable and (if electrically and color neutral) 
frequently makes an excellent DM candidate. Ex- 
amples of WIMP particles come from 1 . lightest 
neutralino state in SUSY theories with conserved 
-R-paritv[il6]. 2. lightest Kaluza-Klein excitations 
from extra-dimensional theories with conserved 
A'if-parity[|l7lll8] and 3. lightest T-odd particles 
in Little Higgs theories with conserved r-parity[ 

muniiiii 11110 

It is possible to calculate the thermal WIMP 
abundance from the Big Bang using very gen- 
eral principles. The initial condition is that at 
early universe temperatures T > mwiMP, the 
WIMPs would have been in thermal equilibrium 
with the cosmic soup. In this case, their abun- 
dance follows straightforwardly from equilibrium 
statistical mechanics. As the universe expands 
and cools, ultimately the WIMPs fall out of ther- 
mal equilibrium at a temperature where the ex- 
pansion rate of the universe equals the WIMP an- 
nihilation rate, because then the WIMPS are un- 
able to find one another to annihilate fast enough: 
this is known as the freeze-out temperature Tp. 
As a result, the WIMP density does not drop ex- 
ponentially as the Universe continues to cool, but 
reduces only as due to the expansion of the 
Universe. The WIMP abundance after freeze-out 
can be found by solving the Boltzmann equation 
in a Friedman-Robertson- Walker universe for the 
WIMP number density. The WIMP mass density 



today, p{Tq), is then given by 



-"^We point out that it has recently been argued [,23, that T-parity 
is generically not conserved because of anomahes in the quantum 
theory. It has, however, been pointed out that whether T-parity 
is or is not conserved can only be definitively addressed only in 
the context of a UV-completion of the model [ 2T. 
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where ~ 2.72 K is the current temperature 
of the CMB, To is the corresponding neutrahno 
temperature, t;, ^ 100 is the number of rela- 
tivistic degrees of freedom at WIMP freeze-out, 
(ctw) is the thermally averaged WIMP annihi- 
lation cross section times relative velocity, and 
xp = 7f/wwimp — 1/20 is the scaled freeze- 
out temperature. But for the fact that photons 
are reheated as various species decouple, the tem- 
peratures of the WIMPs and photons would have 
been the same. Since the reheating process is as- 
sumed to be isentropic, the ratio is sim- 
ply given by the ratio of the number of effec- 
tive degrees of freedom at freeze-out to that to- 
day, and is about 20. Dividing by the closure 
density pc = 8.1 x 10"^"^/!^ GeV"' then gives us 
^wiMph^, where h is the Hubble parameter in 
units of 100 km/s/Mpc. For s-wave annihilation, 
{(jv) is independent of x; then, for fi/i^ ^ 0.1, 
we find it is about 10 pb ~ about the size of an 
electroweak cross section for annihilation of non- 
relativistic particles with a mass of about 50 GeV, 
not far from the weak scale! This provides inde- 
pendent astrophysical evidence that new physics 
- the dark matter particle - may well be lurking 
at the weak scale! The co-incidence of the scale 
of dark matter with the scale of EWSB is some- 
times referred to as the WIMP miracle, and sug- 
gests that the new physics that governs EWSB 
may coincide with the DM sector, and inspires 
many to believe that WIMPs are the prime can- 
didate to constitute the cold dark matter of the 
universe 

• SuperWIMPs: Super WIMPS are electrically and 
color neutral stable DM candidates that inter- 
act with much smaller strength (perhaps only 
gravitationally) than WIMPS. Such particles of- 
ten occur in particle physics theories that include 
WIMPs. Examples include 1. the lightest n = 1 
level KK graviton Gj^^ in extra-dimensional the- 
ories, 2. the gravitino G (the superpartner of 
the graviton) in SUSY theories and 3. the ax- 
ino a (the fermionic member of the axion super- 
multiplet). Since superWIMP interactions with 



2See, however, Ref. [ [25) . 
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ordinary matter have strengths far below conven- 
tional weak interaction strengths, they are not 
expected to yield observable signals in DD or ID 
search experiments. However, they can lead to in- 
triguing new phenomena at collider experiments 
such as LHC and ILC. If every WIMP decays 
to a superWIMP, then superWIMPs inherit the 
thermally produced number density of WIMPs, 
and their contribution to flcDMh^ is reduced 
from the corresponding would-be WIMP contri- 
bution by the ratio of the superWIMP to WIMP 
masses. The superWIMPs produced from WIMP 
decay may be either warm or cold dark matter 
depending on the WIMP lifetime and WIMP- 
superWIMP mass gap [US]. SuperWIMPs may 
also be produced during the re-heating of the Uni- 
verse after inflation; this component of their relic 
abundance is cold, and its magnitude depends on 
the reheating temperature Tr. 

Of the possibilities mentioned above, supersymmetry 
stands out for several reasons. Weak scale supersym- 
metry provides an elegant mechanism to stabilize the 
weak scale against runaway quantum corrections to the 
Higgs scalar mass that arise when the SM is embedded 
into a larger theory that includes particles with masses 
hierarchically larger than the weak scale, e.g. grand 
unified theories (GUTs). Unless the Higgs boson mass 
parameter is tuned with uncanny precision, these cor- 
rections drive the weak scale as well as the physical 
Higgs boson mass to the GUT scale. The supersym- 
metric extension of the SM, with weak scale superpart- 
ners requires no such a fine tuning, and (unlike many 
examples discussed above) provides a framework that 
is perturbatively valid all the way up to the GUT or 
Planck scale. 

SUSY theories thus naturally meld with GUTs, pre- 
serving many of their successes, and providing suc- 
cessful predictions where non-SUSY GUTS appear to 
fail. The latter include the celebrated unification of 
gauge couplings and the value of the ratio m^/m^. 
In many SUSY models with unified values of scalar 
mass parameters renormalized at an ultra-high energy 
scale, radiative corrections drive the weak scale squared 
Higgs boson mass parameter to negative values trigger- 
ing EWSB if the top quark mass is in the range 150- 
200 GeV. This radiative EWSB mechanism was discov- 
ered in the mid-1980s, well before the top mass was 
determined to be '--^ 172 GeV by experiments at the 
Fermilab Tevatron. In addition, fits to precision elec- 
troweak measurements - plotted on the rrit vs. Mw 
plane - now indicate a slight preference for SUSY (with 



light sparticles) over the SM['5T. 

Although weak scale SUSY theories have the very 
attractive features noted above, the presence of many 
new scalar fields also gives rise to potential new prob- 
lems not present in the SM. If supersymmetry is bro- 
ken in an ad hoc manner, flavour-changing processes 
(that do not also change electric charge) occur at un- 
acceptably large rates, as do some CP-violating pro- 
cesses. This is probably a clue about the (presently un- 
known) mechanism by which the superpartners acquire 
SUSY-breaking masses. But the most severe problem 
caused by the appearance of scalars is that we can write 
renormalizable interactions that violate baryon and/or 
lepton number conservation. These interactions would 
cause the proton to decay within a fraction of a second, 
in sharp contrast to a lower limit on its life-time in ex- 
cess of 10^^ years (independent of the mode of decay)! 
To forbid these potentially disastrous interactions, we 
need to posit an additional conservation law, which is 
often taken to be the conservation of a parity-like quan- 
tum number (referred to as i?-parity) taken to be -1-1 
for ordinary particles and —1 for their SUSY partners. 
As a result, the lightest SUSY particle must be stable 
(since all lighter particles have R—+1). 

Unlike the SM, SUSY theories with a conserved R- 
parity naturally include several candidates for DM. All 
that is needed is that the lightest superpartner be elec- 
trically and color neutral. These include, but are not 
limited to: 1. the lightest neutralino Zi, a true WIMP 
candidate, 2. the gravitino G, a gravitationally inter- 
acting spin- 1 superWIMP candidate, 3. the spin-i ax- 
ino a, which is the superpartner of the axion, and 4. the 
superpartner of a sterile neutrino. The super-partner of 
ordinary neutrinos is excluded as galactic DM because 
it would already have been detected by direct searches 
for DM. The axino interaction strength is between that 
of a true WIMP and a gravitino superWIMP. 

Finally, we remark here that the SM does not include 
a viable mechanism for baryogenesis in the early uni- 
verse, primarily because the CP violation is too small. 
In SUSY theories, with their added richness, several 
mechanisms appear to be possible: electroweak baryo- 
genesis, leptogenesis (which is connected to GUT the- 
ories and neutrino mass), so-called Affleck-Dine baryo- 
genesis involving decay of flat directions of the SUSY 
scalar potential and finally, the possibility of infiaton 
decay to heavy neutrino states. 

Despite the lack of direct evidence for SUSY, its 
many attractive features lead many theorists to expect 
weak scale supersymmetry to manifest itself as the next 
paradigm for the laws of physics. While SUSY could 
have fortuitously revealed itself in experiments at LEP 
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or the Tevatron, the LHC is the first facility designed 
to directly probe the weak scale energy regime where 
superpartners are naturally expected. We will, for the 
most part, discuss supersymmetric theories in the re- 
mainder of this article and show that data from the 
LHC as well as from other DD and ID experiments will 
incisively test the weak scale SUSY idea. We will briefly 
return to other ideas with non-SUSY WIMPS in Sec. [71 

4. Supersymmetric theories 

The representations of the SM make a clear dis- 
tinction between the "matter" and "force" sectors of 
the theory. The spin-half matter particles have differ- 
ent gauge quantum numbers from the spin-one gauge 
bosons (which necessarily must be in the adjoint repre- 
sentation of the gauge group) that mediate the strong 
and electroweak interactions. Spin-zero fields, which 
are essential for spontaneous EWSB (and which medi- 
ate a non-gauge force between particles), belong to yet 
another representation. In supersymmetric theories, 
where bosons and fcrmions belong to the same super- 
multiplet, bosons and fermions transform the same way, 
providing a level of synthesis never previously attained. 

The superfield formalism, where bosonic and 
fermionic fields are combined into a single superfield, 
provides a convenient way for constructing supersym- 
metric models of particle physics. This is analogous to 
the familiar isospin formalism where particles of differ- 
ent charge are combined to form an isomultiplet. Chi- 
ral scalar superfields include one chiral-componcnt of a 
spin-half fermion, together with a complex scalar field, 
the superpartner of this chiral fermion. A massive Dirac 
fermion necessarily has two chiral components, and so 
needs two chiral superfields to describe it. For example, 
the Dirac electron therefore has two complex scalar su- 
perpartners (denoted by and en), one corresponding 
to each chirality of the electron/positron. Notice that 
the number of polarization states for fermions (four, 
because there are two polarizations each for the elec- 
tron and positron) is exactly the same as the number 
of bosonic polarization states (each complex spin zero 
field corresponds to two polarization states, one for the 
spin-zero particle, and one for the spin-zero antiparti- 
cle) . This equality of bosonic and fermionic degrees of 
freedom is a general feature of SUSY models. Moreover, 
the gauge quantum numbers for the spin-zero partners 
of the chiral fermion fields must be the same as for 
the corresponding fermions, so that the usual minimal 
coupling prescription completely fixes the gauge inter- 
actions of these particles. 

Gauge superfields include spin-1 gauge bosons along 



with spin-i self-conjugate (or Majorana) gauginos, 
both tranforming under the adjoint representation. Fi- 
nally, there are gravitational supermultiplets containing 
massless spin-2 graviton fields and spin-| gravitinos. 
These are all representations of = 1 supersymme- 
try, where there is just one super-charge. We will focus 
here only on iV = 1 SUSY since it leads most directly to 
phenomenologically viable models with chiral fermions. 

The superfield formalism[ 28, 2£, J50J facilitates the 
construction of a supersymmetric version of the Stan- 
dard Model, known as the Minimal Supersymmetric 
Standard Model, or MSSM. As explained above, for 
each quark and lepton of the SM, the MSSM necessarily 
includes spin-0 superpartners and along with 
and ^fl, whose gauge quantum numbers are fixed to be 
the known gauge quantum numbers of the correspond- 
ing fermions. Thus, for example, the right-handed up 
quark scalar (usually denoted by uji) is a color-triplet, 
weak isosinglet with the same weak hypercharge 4/3 as 
the right-handed up-quark. The MSSM thus includes 
a plethora of new scalar states: e^, cr, v^Li ul, ur, 
d-L, dfi in the first generation, together with analogus 
states for the other two generations. Spin-zero squark 
partners of quarks with large Yukawa couplings undergo 
left-right mixing: thus, the t^ and tji states mix to form 
mass eigenstates - ti and ^2 ^ ordered from lowest to 
highest mass. 

The spin-0 Higgs bosons are embedded in Higgs su- 
perfields, so that the MSSM also includes spin-i higgsi- 
nos. Unlike in the SM, the same Higgs doublet cannot 
give a mass to both up- and down- type fermions with- 
out catastrophically breaking the underlying supersym- 
metry. Thus the MSSM includes two Higgs doublets 
instead of one as in the SM. This gives rise to a richer 
spectrum of physical Higgs particles, including neutral 
light h and heavy H scalars, a pseudoscalar A and a 
pair of charged Higgs bosons H"^ . 

The gauge sector of the MSSM contains gauge bosons 
along with spin-half gauginos in the adjoint represen- 
tation of the gauge group: thus, along with eight col- 
ored gluons, the MSSM contains eight colored spin-i 
gluinos. Upon electroweak symmetry breaking, the four 
gauginos of SU{2)l x U{1)y mix (just as the SU{2)l 
and U{1)y gauge bosons mix) amongst themselves and 
the higgsinos, to form charginos - and ~ and 
neutralinos - Zi, Z2, and Z^. The Zi state, the 
lightest neutralino, is often the lightest supersymmetric 
particle (LSP), and turns out to be an excellent WIMP 
candidate for CDM in the universe. 

If nature is perfectly supersymmetric, then the spin- 
superpartners would have exactly the same mass as 
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the corresponding fermions. Charged spin-0 partners 
of the electron with a mass of 0.51 MeV could not have 
evaded experimental detection. Their non-observation 
leads us to conclude that SUSY must be a broken sym- 
metry. In the MSSM, SUSY is broken explicitly by 
including so-called soft SUSY breaking (SSB) terms in 
the Lagrangian. The SSB terms preserve the desirable 
features of SUSY, such as the stabilization of the scalar 
sector in the presence of radiative corrections, while 
lifting the superpartner masses in accord with what is 
necessary from experiment. It is important to note that 
the equality of dimensionless couplings between par- 
ticles and their superpartners is still preserved (mod- 
ulo small effects of radiative corrections): in particu- 
lar, phenomenologically important gauge interactions 
of superpartners and the corresponding interactions of 
gauginos remain (largely) unaffected by the SSB terms. 

The addition of the SSB Lagrangian terms may seem 
ad-hoc and ugly. It would be elegant if instead super- 
symmetry could be spontaneously broken. But it was 
recognized in the early to mid-1980's that models where 
global SUSY is spontaneously broken at the weak scale 
ran into serious difficulties. The situation is very dif- 
ferent if we elevate SUSY from a global symmetry to 
a local one. In local SUSY, we are forced to include 
the graviton/gravitino super-multiplet into the theory, 
in much the same way that we have to include spin-1 
gauge fields to maintain local gauge invariance of Yang- 
Mills theories. Theories with local SUSY are known as 
supergravity (SUGRA) theories because they are super- 
symmetric and necessarily include gravity. Moreover, 
the gravitational sector of the theory reduces to general 
relativity in the classical limit. Within the framework of 
SUGRA it is possible to add an additional sector whose 
dynamics spontaneously breaks SUSY but which inter- 
acts with SM particles and their superpartners only via 
gravity (the so-called hidden sector). The spontaneous 
breakdown of supersymmetry results in a mass for the 
gravitino in the same way that in local gauge theo- 
ries gauge bosons acquire mass by the Higgs mecha- 
nism. This is, therefore, referred to as the super-Higgs 
mechanism. The remarkable thing is that because of 
the gravitational couplng between the hidden and the 
MSSM sectors, the effects of spontaneous supersymme- 
try breaking in the hidden sector are conveyed to the 
MSSM sector, and (provided the SUSY-breaking scale 
in the hidden sector is appropriately chosen) weak scale 
SSB terms that lift the undesirable degeneracies be- 
tween the masses of SM particles and their superpart- 
ners are automatically induced. Indeed, in the limit 
where Mpi — )■ oo (keeping the gravitino mass fixed) , we 
recover a global SUSY theory along with the desired 



SSB terms! The gravitino typically has a weak scale 
mass and decouples from particle physics experiments 
because of its tiny gravitational couplings. For reasons 
that we cannot discuss here, these locally supersymmet- 
ric models are free[[28l [29j [30] of the above-mentioned 
difficulties that plague globally supersymmetric mod- 
els. 

Motivated by the successful unification of gauge cou- 
phngs at a scale Mgut 2 x 10^'^ GeV in the MSSM, we 
are led to construct a GUT based on local supersym- 
metry. In this case, the theory renormalized at Q = 
Mqjjt contains just one gaugino mass parameter mi/2- 
Renormalization effects then split the physical gaug- 
ino masses in the same way the measured values of the 
gauge couplings arise from a single unified GUT scale 
gauge coupling. In general, supergravity models give 
rise to complicated mass matrices for the scalar super- 
partners of quarks and leptons, with concomitant flavor 
violation beyond acceptable levels. However, in models 
with universal soft SUSY breaking terms, a super-GIM 
mechanism suppresses flavor violating processes [[31]. In 
what has come to be known as the minimal supergravity 
(mSUGRA) model, a universal scalar mass toq a-nd also 
a universal SSB scalar coupling Aq are assumed to exist 
at a high scale Q = A^gut — Afpi- The physical masses 
of squarks and sleptons are split after renormalization, 
and can be calculated using renormalization group tech- 
niques. Typically, in the mSUGRA model, we have 
mq ^ m^^ ^ ^Ir- Although the Higgs scalar mass pa- 
rameters also start off at the common value toq at the 
high scale, the large value of the top quark Yukawa cou- 
pling drives the corresponding squared mass parameter 
to negative values and EWSB is radiatively broken as 
we have already discussed. Within this framework, the 
masses and couplings required for phenomenology are 
fixed by just a handful of parameters which are usually 
taken to be, 

mo, mi/2, Aq, tan/3, and sign{fi). (2) 

Here tan f3 is the ratio of the vacuum expectation values 
of the Higgs fields that give masses to up and down type 
fermions, and fi is the supersymmetric higgsino mass 
parameter whose magnitude is fixed to reproduce the 
measured value of Mz- If all parameters are real, then 
potentially large CP- violating effects are suppressed as 
well. Computers codes such as Isajet, SuSpect, Soft- 
SUSY and Spheno that calculate the full spectrum of 
sparticle and Higgs boson masses are publicly available [ 

The mSUGRA model (sometimes referred to as the 
constrained MSSM or CMSSM) serves as a paradigm 
for many SUSY phenomenological analyses. However, 
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it is important to remember that it is based on many 
assumptions that can be tested in future coUider ex- 
periments but which may prove to be incorrect. For 
instance, in many GUT theories, it is common to get 
non-universal SSB parameters. In addition, there are 
other messenger mechanisms besides gravity. In gauge- 
mediated SUSY breaking models (GMSB)[03], a spe- 
cial messenger sector is included, so gravitinos may be 
much lighter than all other sparticles, with implications 
for both collider physics and cosmology. In anomaly- 
mediated SUSY breaking (AMSB) models['34], gravita- 
tional anomalies induce SSB terms, and the gravitino 
can be much heavier than the weak scale. There are 
yet other models [ 35 where SSB parameters get com- 
parable contributions from gravity-mediated as well as 
from anomaly-mediated sources, and very recently, also 
from gauge-mediation[ I36j . The pattern of superpart- 
ner masses is sensitive to the mediation-mechanism, 
so that we can expect collider experiments to reveal 
which of the various mechanisms that have been pro- 
posed are actually realized in nature. We also mention 
that in both the GMSB and AMSB models, it is some- 
what less natural (but still possible!) to obtain the 
required amount of SUSY dark matter in the Universe. 
Although these are all viable scenarios, they have not 
been as well scrutinized as the mSUGRA model. 

5. Supersymmetric dark matter 

5.1. Neutralino relic density 

Once a SUSY model is specified, then given a set 
of input parameters, it is possible to all compute su- 
perpartner masses and couplings necessary for phe- 
nomenology. We can then use these to calculate scat- 
tering cross sections and sparticle decay patterns to 
evaluate SUSY signals (and corresponding SM back- 
grounds) in collider experiments. We can also check 
whether the model is allowed or excluded by experi- 
mental constraints, either from direct SUSY searches, 
e.g. at LEP2 which requires that > 103.5 GeV, 

ms ~ 100 GeV, and ruh > 114.4 GeV (for a SM-hke 
light SUSY Higgs boson h), or from indirect searches 
through loop effects from SUSY particles in low energy 
measurements such as B{b s-f) or {g — 2)^. We can 
also calculate the expected thermal LSP relic density. 
To begin our discussion, we will first assume that the 
lightest neutralino Zi is the candidate DM particle. 

As mentioned above, the relic density calculation in- 
volves solving the Boltzmann equation, where the neu- 
tralino density changes due to both the expansion of 
the Universe and because of neutralino annihilation 
into SM particles, determined by the thermally aver- 



aged ZiZi annihilation cross section. An added com- 
plication occurs if neutralino co- annihilation is possi- 
ble. Co-annihilation occurs if there is another SUSY 
particle close in mass to the Zi, whose thermal relic 
density (usually suppressed by the Boltzmann factor 
exp ~'^^^ ) is also significant. In the mSUGRA model, 
co-annihilation may occujt^from a stau, fi, a stop ti 
or the lighter chargino Wi. For instance, in some 
mSUGRA parameter-space regions the fi and Zi are 
almost degenerate, so that they both have a signifi- 
cant^ density in the early universe, and reactions such 
as Zifi T"f occur. Since the electrically charged fi 
can also annihilate efficiently via electromagnetic inter- 
actions, this process also alters the equilibrium density 
of neutralinos. All in all, there are well over a thou- 
sand neutralino annihilation and co-annihilation reac- 
tions that need to be computed, involving of order 7000 
Feynman diagrams. There exist several publicly aval- 
able computer codes that compute the neutralino relic 
density: these include DarkSUSY[[32, MicroMegas[l38] 
and IsaReD[l39] (a part of the Isatools package of Isajet[ 

mi 

As an example, we show in Fig. [3] the mg vs. ?ni/2 
plane from the mSUGRA model, where we take = 0, 

> 0, mt = 171.4 GeV and tan/3 = 10. The red- 
shaded regions are not allowed because either the fi 
becomes the lightest SUSY particle, in contradiction 
to negative searches for long lived, charged relics (left 
edge), or EWSB is not correctly obtained (lower-right 
region). The blue-shaded region is excluded by LEP2 
searches for chargino pair production (m^ < 103.5 
GeV). We show contours of squark (solid) and gluino 
(dashed) mass (which are nearly invariant under change 
of Aq and tan/3). Below the magenta contour near 
"^1/2 ^ 200 GeV, mil < 110 GeV, which is roughly 
the LEP2 lower limit on mh in the model. The thin 
green regions at the edge of the unshaded white re- 
gion has fi— : 0.094 — 0.129 where the neutralino 
saturates the observed relic density. In the adjoining 
yellow regions, ft- h? < 0.094, so these regions require 
multiple DM components. The white regions all have 
fi— > 0.129 and so give too much thermal DM: they 
are excluded in the standard Big Bang cosmology. 

The DM-allowed regions are classified as follows: 

• At very low mo and low mi ^2 values is the so- 
called bulk annihilation region[l4T] . Here, sleptons 
are quite light, so ZiZi ££ via t-channel slep- 
ton exchange. In years past (when flcBMh^ 0.3 
was quite consistent with data) , this was regarded 
as the favored region. But today LEP2 sparticle 
search limits have increased the LEP2-forbidden 
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mSUGRA : tan|3=10, A„=0, ^>0, m,=171.4 GeV 




0/2345 



m„ (TeV) 

0<£2h^< 0.094 • Excluded 

• 0.094 < ah" < 0.129 • LEP2 



Figure 3. DM-allowed regions m the mo — mi/2 plane of 
the mSUGRA model for tan P = 10 with Aq = and fj,> 0. 



region from below, while the stringent bound 
^CDMh^ < 0.13 has pushed the DM-allowed re- 
gion down. Now hardly any bulk region survives 
in the mSUGRA model. 

• At low Too and moderate toi/2, there is a thin 
strip of (barely discernable) allowed region adja- 
cent to the stau-LSP region where the neutralino 
and the lighter stau were in thermal equilibrium 
in the early universe. Here co-annihilation with 
the light stau serves to bring the neutralino relic 
density down to its observed value [[42]. 

• At large toq, adjacent to the EWSB excluded re- 
gion on the right, is the hyperbolic branch/focus 
point (HB/FP) region, where the superpotential 

parameter becomes small and the higgsino- 
content of Zi increases significantly. Then Zi 
can annihilate efficiently via gauge coupling to its 
higgsino component and becomes mixed higgsino- 
bino DM. If to- > Mw, Mz, then ZiZi 
WW, ZZ, Zh is enhanced, and one finds the cor- 
rect measured relic density [ 43^. 

We show the corresponding situation for tan (3 — 52 
in Fig. m While the stau co-annihilation and the 
HB/FP regions are clearly visible, we see that now a 
large DM consistent region now appears. 



• In this region, the value of niA is small enough 
so that ZiZi can annihilate into bb pairs through 
s-channel A (and also H) resonance. This re- 
gion has been dubbed the A- funnel [ [44]. It can 
be quite broad at large tan/3 because the width 
Ta can be quite wide due to the very large b- 
and T- Yukawa couplings. If tan/3 is increased 
further, then ZiZi annihilation through the (vir- 
tual) A* is large all over parameter space, and 
most of the theoretically-allowed parameter space 
becomes DM-consisten. For even higher tan /3 val- 
ues, the parameter space collapses due to a lack 
of appropriate EWSB. 



mSUGRA : A„ = 0, |i > 0, tan|3 = 52, m, = 171.4 GeV 
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0<Qh'< 0.094 • LEP2 Excl. : m^, < 103.5 GcV 

< 0.094 <£2h-< 0.129 — m|,= llOGeV 

Figure 4. DM-allowed regions tn the mo — mi/2 plane of 
the mSUGRA model for tan/3 = 52 wtth Aq = and fi > 0. 
The various colors of shading is as in Fig. 



It is also possible at low mi/2 values that a light 
Higgs h resonance annihilation region can occur just 
above the LEP2 excluded region[l45]. Finally, if is 
large and negative, then the ti can become light, and 
rrif^ TO— , so that stop-neutralino co- annihilation [ H5] 
can occur. 

Up to now, we have confined our discussion to the 
mSUGRA framework in which compatibility with ([T]) is 
obtained only over selected portions of the toq — tox/2 
plane. The reader may well wonder what happens if we 
relax the untested universality assumptions that un- 
derlie mSUGRA. Without going into details, we only 
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mention here that in many simple one-parameter ex- 
tensions of mSUGRA where the universahty of mass 
parameters is relaxed in any one of the matter scalar, 
the Higgs scalar, or the gaugino sectors, all points in 
the mg — mi/2 plane become compatible with the relic 
density constraint due to a variety of mechanisms: these 
are catalogued in Ref. [ [3H] ■ Implications of the relic 
density measurement for collider searches must thus be 
drawn with care. 

5.2. Neutralino direct detection 

Fits to galactic rotation curves imply a local relic 
density of pcoM ^ 0.3 GeV/cm^. For a 100 GeV 
WIMP, this translates to about one WIMP per coffee 
mug volume at our location in the galaxy. The goal 
of DD experiments is to detect the very rare WIMP- 
nucleus collisions that should be occuring as the earth, 
together with the WIMP detector, moves through the 
DM halo. 

DD experiments are usually located deep under- 
ground to shield the experimental apparatus from back- 
ground due to cosmic rays and ambient radiation from 
the environment or from radioactivity induced by cos- 
mic ray exposure. One technique is to use cryogenic 
crystals cooled to near absolute zero, and look for 
phonon and ionization signals from nuclei recoiling from 
a WIMP collision. In the case of the CDMS experiment [ 
I49j at the Soudan iron mine, target materials include 
germanium and silicon. Another technique uses noble 
gases cooled to a liquid state as the target. Here, the 
signal is scintillation light picked up by photomultiplier 
tubes and ionization. Target materials include xenon[ 
150] . argon and perhaps neon. These noble liquid de- 
tectors can be scaled up to large volumes at relatively 
low cost. They have the advantage of fiducialization, 
wherein the outer layers of the detector act as an ac- 
tive veto against cosmic rays or neutrons coming from 
phototubes or detector walls: only single scatters from 
the inner fiducial volume qualify as signal events. A 
third technique, typified by the COUPP experiment [ 
IST] . involves use of superheated liquids such as CF^I 
located in a transparent vessel. The nuclear recoil from 
a WIMP-nuclcon collision then serves as a nucleation 
site, so that a bubble forms. The vessel is monitored 
visually by cameras. Background events are typically 
located close to the vessel wall, while neutron interac- 
tions are likely to cause several bubbles to form, instead 
of just one, as in a WIMP collision. This technique al- 
lows for the use of various target liquids, including those 
containing elements such as fluorine, which is sensitive 
to spin- dependent interactions. 

The cross section for WIMP-nucleon collisions can be 



calculated, and in the low velocity limit separates into 
a coherent spin-independent component (from scatter- 
ing mediated by scalar quarks and scalar Higgs bosons) 
which scales as nuclear mass squared, and a spin- 
dependent component from scattering mediated by the 
Z boson or by squarks, which depends on the WIMP 
and nuclear spins[ 29 . The scattering cross section per 
nucleon versus tovf/mp serves as a figure of merit and 
facilitates the comparison of the sensitivity of various 
experiments using different target materials. 



Spin-independent Direct Detection 
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Figure 5. The spin-independent neutralino-proton scatter- 
ing cross-section vs m~ in a variety of SUSY models, com- 
patible with collider constraints where thermally produced 
Big Bang neutralinos saturate the observed dark matter den- 
sity. 



In Fig. [51 we show the spin-independent Zip cross 

section versus m~ for a large number of SUSY mod- 

^1 

els (including mSUGRA) . Every color represents a dif- 
ferent model. For each model, parameters are chosen 
so that current collider constraints on sparticle masses 
are satisfied, and further, that the lightest neutralino 
(assumed to be the LSP) saturates the observed relic 
abundance of CDM. Also shown is the sensitivity of 
current experiments together with projected sensitivity 
of proposed searches at superCDMS, Xenon-100, LUX, 
WARP and at a ton-sized noble liquid detector. The 
details of the various models are unimportant for our 
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present purpose. The key thing to note is that while 

the various models have a branch where ctsi(p-^i) islls 

off with m~ , there is another branch where this cross- 
^1 

section asymptotes to just under 10^® Db[ ITfl H51 [52] , 
Points in this branch (which includes the HB /FP region 
of mSUGRA), are consistent with ([T]) because Zi has 
a significant higgsino component. Neutralinos with an 
enhanced higgsino content can annihilate efffciently in 
the early universe via gauge interactions. Moreover, 
since the spin-independent DD amplitude is mostly 
determined by the Higgs boson-higgsino-gaugino cou- 
pling, it is large in models with MHDM which has 
both gaugino and higgsino components. Thus the en- 
hanced higgsino component of MHDM increases both 
the neutralino annihilation in the early universe as well 
as the spin-independent DD rate. The exciting thing is 
that the experiments currently being deployed- such as 
Xenon-100, LUX and WARP- will have the sensitivity 
to probe this class of models. To go further will require 
ton-size or greater target material. 

We note here that if towimp ~ 150 GeV, then it 
may be possible to extract the WIMP mass by mea- 
suring the energy spectrum of the recoiling nuclear tar- 
gets[|54l. Typically, of order 100 or more events are 
needed for such a determination to 10-20%. For higher 
WIMP masses, the recoil energy spectrum varies lit- 
tle, and WIMP mass extraction is much more difficult. 
Since the energy tranfer from the WIMP to a nucleus 
is maximized when the two have the same mass, DD 
experiments with several target nuclei ranging over a 
wide range of masses would facilitate the distinction 
between somewhat light and relatively heavy WIMPs, 
and so, potentially serve to establish the existence of 
multiple WIMP components in our halo. 

5.3. Indirect detection of neutralinos 

As explained in Sec. [U there are also a number of in- 
direct WIMP search techniques that attempt to detect 
the decay products from WIMP annihilation at either 
the center of the sun, at the galactic center, or within 
the galactic halo. 

5.3.1. Neutrino telescopes 

Neutrino telescopes such as ANTARES or IceCube 
can search for high energy neutrinos produced from 
WIMP- WIMP annihilation into SM particles in the 
core of the sun (or possibly the earth). The technique 
involves detection of multi-tens of GeV muons produced 
by interactions with polar ice (IceCube) or ocean wa- 
ter (ANTARES). The muons travel at a speeds greater 
than the speed of light in the medium, thus leaving a 
tell-tale signal of Cerenkov light which is picked up by 



arrays of phototubes. The IceCube experiment, cur- 
rently being deployed at the south pole, will monitor a 
cubic kilometer of ice in search oi ^ fi conversions. 
It should be fully deployed by 2011. The experiment is 
mainly sensitive to muons with E"^ > 50 GeV. 

In the case of neutralinos of SUSY, mixed higgsino 
dark matter (MHDM) has a large (spin-dependent) 
cross-section to scatter from hydrogen nuclei via Z- 
exchange and so is readily captured. Thus, in the 
HB/FP region of mSUGRA, or in other SUSY models 
with MHDM, we expect observable levels of signal ex- 
ceeding 40 events/km^ /yr with > 50 GeV. For the 
mSUGRA model, the IceCube signal region is shown 
beneath the magenta contour labelled /i in Fig. [6j'55'. 
These results were obtained using the Isajet-DarkSUSY 
interface [ [37]. Notice that DD signals are also observ- 
able in much the same region (below the contour la- 
belled DD) where the neutrahno is MHDM. 
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Figure 6. The projected reach of various colliders, di- 
rect and indirect dark matter search experiments in the 
mSUGRA model. For the indirect search results we have 
adopted the conservative default DarkSUSY isotropic DM 
halo density distribution. Plot is from Ref. [\55f. 
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5.3.2. Anti-matter from WIMP halo annihila- 
tions 

WIMP annihilation in the galactic halo offers a differ- 
ent possibility for indirect DM searches. Halo WIMPs 
annihilate equally to matter and anti-matter, so the 
rare presence of high energy anti-matter in cosmic ray 
events - positrons e^, anti-protons p, or even anti- 
deuterons D - offer possible signatures. Positrons pro- 
duced in WIMP annihilations must originate relatively 
close by, or else they will find cosmic electrons to anni- 
hilate against, or lose energy via bremsstrahlung. Anti- 
protons and anti-deuterons could originate further from 
us because, being heavier, they are deflected less and so 
lose much less energy. The expected signal rate depends 
on the WIMP annihilation rate into anti-matter, the 
model for the propogation of the anti-matter from its 
point of origin to the earth, and finally on the assumed 
profile of the dark matter in the galactic halo. Several 
possible halo density profiles are shown in Fig. [T] We 
see that while the local WIMP density is inferred to a 
factor of ^ 2-3 (we are at about 8 kpc from the Galactic 
center), the DM density at the galactic center is highly 
model-dependent close to the core. Since the ID signal 
should scale as the square of the WIMP density at the 
source, positron signals will be uncertain by a factor of 
a few with somewhat larger uncertainty for p and D 
signals that originate further away. Anti-particle prop- 
agation through the not so well known magnetic field 
leads to an additional uncertainty in the predictions. 
The recently launched Pamela space-based anti-matter 
telescope can look for e+ or p events while the balloon- 
borne GAPS experiment will be designed to search for 
anti-deuterons. Anti-matter signals tend to be largest 
in the case of SUSY models with MHDM or when neu- 
tralinos annihilate through the A-resonance [ (56] . 

5.3.3. Gamma rays from WIMP halo annihila- 
tions 

As mentioned in the Introduction, high energy 
gamma rays from WIMP annihilation offer some ad- 
vantages over the signal from charged antiparticles. 
Gamma rays would point to the source, and would 
degrade much less in energy during their journey to 
us. ^ This offers the possibility of the line signal from 
ZiZi —^ 77 processes that occur via box an triangle di- 
agrams. While this reaction is loop-suppressed, it yields 
monoenergetic photons with ~ towimp, and so can 
provide a measure of the WIMP mass. Another possi- 
bility is to look for continuum gamma rays from WIMP 
annihilation to hadrons where, for instance, the gamma 
is the the result of vr" decays. Since the halo WIMPS 
are essentially at rest, we expect a diffuse spectrum of 
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Figure 7. Various predictions for the DM halo in the Milky 
Way as a function of distance from the galactic center. The 
earth is located at r ~ 8 kpc. 



gamma rays, but with < mwiMP- Because gamma 
rays can traverse large distances, a good place to look 
at is the galactic center, where the WIMP density (see 
Fig. [7]) is expected to be very high. Unfortunately, the 
density at the core is also very uncertain, making pre- 
dictions for the gamma ray fiux uncertain by as much as 
four orders of magnitude. Indeed, detection of WIMP 
halo signals may serve to provide information about the 
DM distribution in our galaxy. 

Anomalies have been reported in the cosmic gamma 
ray spectrum. In one example, the Egret experiment [ 
[57] sees an excess of gamma rays with E.^ > 1 GeV. 
Explanations for the Egret GeV anomaly range from 
ZiZi ^ 66 ^ 7 with m~ ~ 60 GeV[[SH], to mis- 
calibration of the Egret calorimeter [[55]. The GLAST 
gamma ray observatory is scheduled for lift-off in 2008 
and should help resolve this issue, as will the upcoming 
LHC searches ['60^,. 

5.4. Gravitino dark matter 

In gravity-mediated SUSY breaking models, graviti- 
nos typically have weak scale masses and, because they 
only have tiny gravitational couplings, are usually as- 
sumed to be irrelevant for particle physics phenomenol- 
ogy. Cosmological considerations, however, lead to the 
gravitino problem, wherein overproduction of graviti- 
nos, followed by their late decays into SM particles, 
can disrupt the successful predictions of Big Bang nu- 
cleosynthesis. The gravitino problem can be overcome 
by choosing an appropriate range for and a low 
enough re-heat temperature for the universe after in- 
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flation[[CT[ as illustrated in Fig. [51 or by hypothesizing 
that the G is in fact the stable LSP, and thus consti- 
tutes the DM [Eg. 
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Figure 8. An illustration of constraints from Big Bang 
nucleosynthesis which require Tr to be below the various 
curves, for the HB/FP region of the mSUGRA model with 
mo = 2397 GeV, mi/a = 300 GeV, Ao = and tan/3 = 30, 
from Kohri et al.[ 611 to which we refer the reader for more 
details. 



Here, we consider the consequences of a gravitino 
LSP in SUGRA models. If gravitinos are produced 
in the pre-inflation epoch, then their number density 
will be diluted away during inflation. After the uni- 
verse inflates, it enters a re-heating period wherein all 
particles can be thermally produced. However, the cou- 
plings of the gravitino are so weak that though graviti- 
nos can be produced by the particles that do partake of 
thermal equilibrium, gravitinos themselves never attain 
thermal equilibrium: indeed their density is so low that 
gravitino annihilation processes can be neglected in the 
calculation of their relic density. The thermal produc- 
tion (TP) of gravitinos in the early universe has been 
calculated, and including EW contributions, is given by 
the approximate expression (valid for ttiq <^ M, [l63j): 
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where Tr is the re-heat temperature. 

Gravitinos can also be produced by decay of the 
next-to-lightest SUSY particle, the NLSP. In the case 



of a long-lived neutralino NLSP, the neutralinos will 
be produced as usual with a thermal relic abundance 
in the early universe. Later, they will each decay as 
Zi 7(5, ZG or hG. The total relic abundance is 
then 



(4) 



The G from NLSP decay may constitute warm/hot 
dark matter depending in the Zi — G mass gap, while 
the thermally produced G wiU be CDM[f2B]. 

The lifetime for neutralino decay to the photon and 
a gravitino is given by [ |64] , 

487rAf2^2 
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the wino and bino components of the Zi [HH], is the 
reduced Planck mass, and r = m^jm- . Similar for- 
mulae (with different mixing angle and r-dependence) 
hold for decays to the gravitino plus a Z or h boson. 
We see that - except when the gravitino is very much 
lighter than the neutralino as may be the case in GMSB 
models with a low SUSY breaking scale - the NLSP 
decays well after Big Bang nucleosynthesis. Such de- 
cays would inject high energy gammas and/or hadrons 
into the cosmic soup post-nucleosynthesis, which could 
break up the nuclei, thus conflicting with the successful 
BBN predictions of Big Bang cosmology. For this rea- 
son, gravitino LSP scenarios usually favor a stau NLSP, 
since the BBN constraints in this case are much weaker. 

Finally, we remark here upon the interesting inter- 
play of baryogenesis via leptogenesis with the nature of 
the LSP and NLSP. For successful thermal leptogenesis 
to take place, it is found that the reheat temperature 
of the universe must exceed ^ 10^" GeV[[nS]. If this is 
so, then gravitinos would be produced thermally with a 
huge abundance, and then decay late, destroying BBN 
predictions. For this reason, some adherents of leptoge- 
nesis tend to favor scenarios with a gravitino LSP, but 
with a stau NLSP[l66]. 

5.5. Axino dark matter 

If we adopt the MSSM as the effective theory below 
-^GUT, and then seek to solve the strong GP problem 
via the Peccei-Quinn solution [ '13], we must introduce 
not only an axion but also a spin-i axino a into the 
theory. The axino mass is found to be in the range 
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of keV-GeV['BT, but its coupling is suppressed by the 
Peccei-Quinn breaking scale fa, which is usually taken 
to be of order 10^ — 10^^ GeV: thus, the axino interacts 
more weakly than a WIMP, but not as weakly as a 
gravitino. The axino can be an compelling choice for 
DM in the universe [[55], 

Like the gravitino, the axino will likely not be in ther- 
mal equilibrium in the early universe, but can still be 
produced thermally via particle scattering. The ther- 
mal production abundance is given bv[ l68l 169] 

oTP,2 . . 61 /1-108\ /lO" GeVy 



VlOO MeV/ \10^ GeV / ' ^ ' 

where fa is the PQ scale, is a model-dependent color 
anomaly factor that enters only as fa/N, and gs is the 
strong coupHng at the reheating scale. 

Also like the gravitino, the axino can be produced 
non-thermally by NLSP decays, where the NLSP abun- 
dance is given by the standard relic density calculation. 
Thus, 

n,h' = nTPh' + -Jl^n^^sph^ (7) 

rriNLSP 

In this case, the thermally produced axinos will be 
CDM for ma ^ 0.1 MeV[[68], while the axinos pro- 
duced in NLSP decay will constitute hot /warm DM[ 
[55]. Since the PQ scale is considerably lower than the 
Planck scale, the lifetime for decays such as Zi ja 
are of order ~ 0.03 sec- well before BBN. Thus, the ax- 
ino DM scenario is much less constrained than gravitino 
DM. 

Note also that if axinos are the CDM of the universe, 
then models with very large fi— ft,^ ~ 100 — 1000 can be 
readily accommodated, since there is a huge reduction 
in relic density upon Zi decay to the axino. This pos- 
sibility occurs in models with multi-TeV scalars (and 
hence a multi-TeV gravitino) and a bino-like Zi. In 
this case with very large there is no gravitino prob- 
lem as long as the re-heat temperature Tr ~ 10^ — 10* 
GeV. This range of Tr is also what is needed to obtain 
successful non-thermal leptongenesis (involving heavy 
neutrino N production via inflaton decay) [ [70] along 
with the correct abundance of axino dark matter [ [71] . 
A scenario along these lines has been proposed [[72] to 
reconcile Yukawa-unified SUSY models, which usually 
predict a vast over-abundance of neutralino DM, with 
the measured relic density. 



6. SUSY DM at the LHC 

6.1. Sparticle production at the LHC 

Direct production of neutralino dark matter at the 
LHC {pp Z\Z\X, where X stands for assorted 
hadronic debris) is of little interest since the high pT 
final state particles all escape the detector, and there is 
little if anything to trigger an event record. Detectable 
events come from the production of the heavier super- 
partners, which in turn decay via a multi-step cascade 
which ends in the stable LSP. 

In many models, the strongly interacting squarks 
and/or gluinos are among the heaviest states. Unless 
these are extremely heavy, these will have large pro- 
duction cross sections at the LHC. Strong interaction 
production mechanisms for their production include, 1. 
gluino pair production gg, 2. squark pair production qq 
and 3. squark-gluino associated production qg. Note 
here that the reactions involving squarks include a huge 
number of subprocess reactions to cover the many fla- 
vors, types (left- and right-), and also the anti-squarks. 
The various possibilities each have different angular de- 
pendence in the production cross sections [[73], and the 
different flavors/ types of squarks each have different de- 
cay modes [TT. These all have to be kept track of in 
order to obtain a reliable picture of the implications 
of SUSY in the LHC detector environment. Squarks 
and gluinos can also be produced in association with 
charginos and neutralinos[ 75 . Associated gluino pro- 
duction occurs via squark exchange in the t or u chan- 
nels and is suppressed if squarks are very heavy. 

If colored sparticles are very heavy, then electroweak 
production of charginos and neutralinos may be the 
dominant sparticle production mechanism at the LHC. 
The most important processes are pair production of 
charginos, W^Wj where = 1,2, and chargino- 

neutralino production, W^Zj, with i ~ 1,2 and j = 
1 — 4. In models jwitjh^unified GUT_scale gaugino masses 
and large |^|, ZW\W\ and WZiW\ couphngs are large 
so that W\W\ and W1Z2 production occurs at signif- 
icant rates. The latter process can lead to the gold- 
plated trilepton signature at the LHC[ 76J. Neutralino 
pair production {pp —> ZiZjX where i,j — 1 — 4) is 
also possible. This reaction occurs at low rates at the 
LHC unless |^| ~ Mi, 2 (as in the case of MHDM). Fi- 
nally, we mention slepton pair production: £^£^, i)f/ 
and vivi, which can give detectable dilepton signals if 
m^-~ 300 GeV[[77]. 

In Fig. [S] we show various sparticle production cross 
sections at the LHC as a function of mg . Strong inter- 
action production mechanisms dominate at low mass. 
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while electroweak processes dominate at high mass. 
The associated production mechanisms are never dom- 
inant. The expected LHC integrated luminosity in the 
first year of running is expected to be around 0.1 fb^^, 
while several tens of fb~^ of data is expected to be 
recorded in the first several years of operation. The 
ultimate goal is to accumulate around 500-1000 fb^^, 
correponding to 10^ — 10^ SUSY events for ~ 1 TeV. 
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Figure 9. Cross sections for production of various sparticles 
at the LHC. Caugino mass unification is assumed. 



6.2. Sparticle cascade decays 

In i?-parity conserving models, sparticles decay to 
lighter sparticles until the decay terminates in the LSP [ 
Frequently, the direct decay to the LSP is ei- 
ther forbidden or occurs with only a small branching 
fraction. Since gravitational interactions are negligible, 
gluinos can only decay via g — > qq, where the q and q 
can be of any flavor or type. If two body decay modes 
are closed, the squark will be virtual, and the gluino 



will decay via three body modes g — s- qqZi, qq'Wj. If 
squarks are degenerate, and Yukawa coupling effects 
negligible, three-body decays to the wino-like chargino 
and neutralino usually have larger branching fractions 
on account of the larger gauge coupling. If < M2, 
gluinos and squarks may thus decay most of the time 
to the heavier charginos and neutralinos, resulting in 
lengthy cascade decay chains at the LHC. 

Squarks decay always to two-body modes — > qg 
if it is kinematically allowed, or q^ q'Wi, qZj, 
while qji — *■ qZj only, since right-squarks do not cou- 
ple to charginos. Sleptons do not have strong inter- 
actions so cannot decay to gluinos. Their electroweak 
decays are similar to corresponding decays of squarks 
II i'W„ £Zj while £r iZj only. ^ 

Charginos may jdecay via two-body modes: Wi 
WZj, ivi, £i>i, ZWj or even to 4>Wj or H~Zj, where 
(p = h,H, A. If two-body modes are inaccessible, then 
three-body decays dominate: Wi ^jff, where / 
and /' are SM fermions which couple to the W. Fre- 
quently, the decay amplitude is dominatedjDy the vir- 
tual W so that the three-body decays of Wi have the 
same branching fractions as^those of the W. Neu- 
tralinos decay via Zi — > WWj, H^Wj, ZZj, (pZj or 
//. If two body neutralino decays are closed, then 
Zi Zjff, where / are the SM fermions. In some 
models, the branching fraction for radiative decays 
Zi — > Zjj (that only occurs at the one- loop level) may 
be significant [ 175]. The cascade decay modes of neu- 
tralinos depend sensitively on model parameters [[75]. 

If tan /3 is large, then b and r Yukawa coupling effects 
become important, enhancing three body decays of ^, 
Wi and Zj to third generation fermions [[80]. For very 
large values of tan/? these decays can even dominate, 
resulting in large rates for &-jet and r-jet production in 
SUSY events [i81J- 

Finally, the various Higgs bosons can be produced 
both directly and via sparticle cascades at the LHC[ 
[H2] ■ Indeed, it may be possible that h is first discovered 
in SUSY events because in a sample of events enriched 
for SUSY, it is possible to identify h via its dominant 
h ^ bb decays rather than via its sub-dominant decay 
modes, as required for conventional searches [[82]. The 
heavier Higgs bosons decay to a variety of SM modes, 
but also to SUSY particles if these latter decays are 
kinematically allowed, leading to novel signatures such 
as H, Z2Z2 ^4:e + Ef'^^l^. 

The cascade decays terminate in the LSP. In the case 
of a Zi LSP, the Zi is a DM candidate, and leaves its 
imprint via E^^^^. In the case of a weak scale G or a 
LSP, then Zi will decay as discussed above. In these 
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cases, the Zi lifetime is long enough that it decays out- 
side the detector, so one still expects large i?™'"" in 
the collider events. An exception arises for the case of 
super-light gravitinos (with masses in the eV to keV 
range) that are possible in GMSB models: see (O- 
Then, the decay may take place inside inside the detec- 
tor, possibly with a large vertex separation. It is also 
possible that the NLSP is charged and quasi-stable, in 
which case collider events may include highly ionizing 
tracks instead of, or in addition to, E^^^^. 

The decay branching fractions depend on the entire 
spectrum of SUSY particle masses and their mixings. 
They are pre-programmed in several codes: Isaiet[l40]. 
SDECAY[[H1] and SphenoifSS]. 

6.3. Event generation for LHC 

Once sparticle production cross sections and decay 
branching fractions have been computed, it is useful to 
embed these into event generator programs to simulate 
what SUSY collider events will look like at LHC. There 
are several steps involved: 

• Calculate all sparticle pair production cross sec- 
tions. Once all initial and final states are ac- 
counted for, this involves over a thousand indi- 
vidual subprocess reactions. In event generation, 
a particular reaction is selected on a probabilistic 
basis, with a weight proportional to its differential 
cross-section. 

• Sparticle decays are selected probabilistically into 
all the allowed modes in proportion to the corre- 
sponding branching fractions. 

• Initial and final state quark and gluon radiation 
are usually dealt with using the parton shower 
(PS) algorithm, which allows for probabilistic 
parton emission based on approximate collinear 
QCD emission matrix elements, but exact kine- 
matics. The PS is also applied at each step of the 
cascade decays, which may lead to additional jet 
production in SUSY collider events. 

• A hadronization algorithm provides a model for 
turning various quarks and gluons into mesons 
and baryons. Unstable hadrons must be further 
decayed. 

• The beam remnants - proton constituents not 
taking part in the hard scattering - must be show- 
ered and hadronized, usually with an independent 
algorithm, so that energy deposition in the for- 
ward detector region may be reliably calculated. 



At this stage, the output of an event generator pro- 
gram is a listing of particle types and their associated 
four- vectors. The resulting event can then be interfaced 
with detector simulation programs to model what the 
actual events containing DM will look like in the envi- 
ronment of a collider detector. 

Several programs are available, including Isaiet[[iD]. 
Pythia['86] and IIerwig[ '87]. Other programs such as 
Madevent[_88j, CompHEP/CalcHEP[^ and Whizard[ 
[90] can generate various 2 n processes including 
SUSY particles. The output of these programs may 
then be used as input to Pythia or Herwig for show- 
ering and hadronization. Likewise, parton level Isajet 
SUSY production followed by cascade decays can be in- 
put to Pythia and Herwig via the Les Houches Event 
format [ [ST]. 

6.4. Signatures for sparticle production 

Unless colored sparticles are very heavy, the SUSY 
events at the LHC mainly result in gluino and squark 
production, followed by their possibly lengthy cascade 
decays. These events, therefore, typically contain very 
hard jets (from the primary decay of the squark and/or 
gluino) together with other jets and isolated electrons, 
muons and taus (identified as narrow one- and three- 
prong jets), and sometimes also photons, from the de- 
cays of secondary charginos and neutralinos, along with 
E^^^^ that arises from the escaping dark matter par- 
ticles (as well as from neutrinos). In models with a 
superlight gravitino, there may also be additional iso- 
lated photons, leptons or jets from the decay of the 
NLSP. The relative rates for various rt-jet -I- m-lepton 
-I- fc-photon +E^^^^ event topologies is sensitive to the 
model as well as to the parameter values, and so provide 
a useful handle for phenomenological analyses. 

Within the SM, the physics background to the classic 
jets + E^^^^ signal comes from neutrinos escaping the 
detector. Thus, the dominant SM backgrounds come 
from W + jets and Z + jets production, tt produc- 
tion, QCD multijet production (including bb and cc pro- 
duction), WW, WZ, ZZ production plus a variety of 
2 — > n processes which are not usually included in event 
generators. These latter would include processes such 
as titi, ttbb, ttW, WWW, WWZ production, etc. De- 
cays of electroweak gauge bosons and the t-quark are 
the main source of isolated leptons in the SM. Various 
additional effects- uninstrumented regions, energy mis- 
measurement, cosmic rays, beam-gas events- can also 
lead to Eif'^'^ events. 

In contrast to the SM, SUSY events naturally tend to 
have large jet multiplicities and frequently an observ- 
able rate for high multiplicity lepton events with large 
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jTjmiss^ Thus, if one plots signal and background versus 
multiplicity of any of these quantities, as one steps out 
to large multiplicity, the expected SUSY events should 
increase in importance, and even dominate the high 
multiplicity channels in some cases. This is especially 
true of isolated multi-lcpton signatures, and in fact it is 
convenient to classify SUSY signal according to lepton 
multiplicity [El]: 

• zero lepton +jets + E™^^^ events, 

• one lepton +jets + E!^^^^ events, 

• two opposite sign leptons +jets + E^^^^ events 
(OS), 

- same-flavor (OSSF), 

- different flavor (OSDF), 

• two same sign leptons +jets + E^^^^ events (SS), 

• three leptons +jets + E^^^^ events (3£), 

• four (or more) leptons +jets + E^^^^ events {M). 

6.5. LHC reach for SUSY 

Event generators, together with detector simulation 
programs can be used to project the SUSY discovery 
reach of the LHC. Given a specific model, one may first 
generate a grid of points that samples the parameter 
(sub)space where signals rates are expected to vary sig- 
nificantly. A large number of SUSY collider events can 
then be generated at every point on the grid along with 
the various SM backgrounds to the SUSY signal men- 
tioned above. Next, these signal and background events 
are passed through a detector simulation program and 
a jet-finding algorithm is implemented to determine the 
number of jets per event above some Exijet) threshold 
(usually taken to be Erijet) > 50-100 GeV for LHC). 
Finally, analysis cuts are imposed which are designed to 
reject mainly SM BG while retaining the signal. These 
cuts may include both topological and kinematic selec- 
tion criteria. For observability with an assumed inte- 
grated luminosity, we require that the signal exceed the 
chance 5 standard deviation upward fluctuation of the 
background, together with a minimum value of (^ 25%) 
the signal to background ratio, to allow for the fact that 
the background is not perfectly known. For lower spar- 
ticle masses, softer kinematic cuts are used, but for high 
sparticle masses, the lower cross sections but higher en- 
ergy release demand hard cuts to optimize signal over 
background. 

In Fig. \Wi we illustrate the SUSY reach of the LHC 
within the mSUGRA model assuming an integrated 



luminosity of 100 fb~^. We show the result in the 
mo — plane, taking Aq = 0, tan (3 = 10 and /i > 0. 
The signal is observable over background in the corre- 
sponding topology below the corresponding curve. We 
note the following. 

1. Unless sparticles are very heavy, there is an ob- 
servable signal in several different event topolo- 
gies. This will help add confidence that one is 
actually seeing new physics, and may help to sort 
out the production and decay mechanisms. 

2. The reach at low mo extends to mi/2 ~ 
1400 GeV. This corresponds to a reach for rriq ~ 
mg - 3.1 TeV. 

3. At large mo, squarks and sleptons are in the 4 — 5 
TeV range, and are too heavy to be produced at 
significant rates at LHC. Here, the reach comes 
mainly from just gluino pair production. In this 
range, the LHC reach is up to mi/2 ^ 700 GeV, 
corresponding to a reach in nig of about 1.8 TeV, 
and may be extended by ^ 15-20% by 6-jet tag- 
ging [[^3]. 
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Figure 10. The 100 fl-^ fb reach of LHC for SUSY m 
the mSUGRA model. For each event topology, the signal is 
observable below the corresponding contour. 



In Fig. [S] we can see a comparison of the LHC reach 
(notice that it is insensitive to tan/3 and sign(/i)) with 
that of the Tevatron (for clean 3£ events with 10 fb~^), 
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and the proposed e'^e~ International Linear Collider 
(ILC), with y/s = 0.5 or 1 TeV along with various dark 
matter DD and ID search experiments. We remark 
that: 

• While LHC can cover most of the relic density al- 
lowed region, the HB/FP region emerges far be- 
yond the LHC reach. 

• As already noted, the DD and ID experiments 
have the greatest sensitivity in the HB/FP region 
where the neutralino is MHDM. In this sense, DD 
and ID experiments complement LHC searches for 
SUSY. 

• The ILC reach is everywhere lower than LHC, ex- 
cept in the HB/FP region. In this region, while 
gluinos and squarks can be extremely heavy, the 
^ parameter is small, leading to a relatively light 
spectrum of charginos and neutralinos. These 
are not detectable at the LHC because the vis- 
ible decay products are too soft. However, since 
chargino pair production is detectable at ILC 
even if the energy release in chargino decays is 
small, the ILC reach extends beyond LHC in this 
region [[M]. 

Finally, we note here that while the results presented 
above are for the LHC reach in the mSUGRA model, 
the LHC reach (measured in terms of and rriq) tends 
to be relatively insensitive to the details of the model 
chosen, as long as gluino and squark production fol- 
lowed by cascade decays to the DM particle occur. 

6.6. Early discovery of SUSY at LHC without 

Recently, it has been pointed out that a SUSY search 
using the traditional jets + E^^^^ signature may not be 
possible for a while after start-up due to various de- 
tector calibration issues. In this case, it is possible to 
abandon using the E^^^^ cut, and instead require a high 
multiplicity of isolated leptons: SS, OSSF, OSDF, 3i. 
The high lepton multiplicity requirement severely re- 
duces SM background while maintaining large enough 
signal rates. In Ref. [[95], it is claimed an LHC reach 
of rrig ~ 750 GeV is possible with just 0.1 fb~^ of inte- 
grated luminosity, without using an E!^^^^ cut. 

6.7. Determination of sparticle properties 

Once a putative signal for new physics emerges at 
LHC, the next step is to establish its origin. This 
will entail detailed measurements of cross sections and 
distributions in various event topologies to gain in- 
sight into the identity of the new particles being pro- 
duced, their masses, decay patterns, spins, couplings 



(gauge quantum numbers) and ultimately mixing an- 
gles. These measurements are not straightforward 
in the LHC environment because of numerous possi- 
ble SUSY production reactions occurring simultane- 
ously, a plethora of sparticle cascade decay possibilities, 
hadronic debris from initial state radiation and lack of 
invariant mass reconstruction due to the presence of 
E^^^^. All these lead to ambiguities and combinatoric 
problems in reconstructing exactly what sort of signal 
reactions are taking place. In contrast, at the ILC, the 
inital state is simple, the beam energy is tunable and 
beam polarization can be used to select out specific 
processes. 

While it seems clear that the ILC is better suited for 
a systematic program of precision sparticle measure- 
ments, studies have shown (albeit in special cases) that 
interesting measurements are also possible at the LHC. 
We go into just a subset of all details here in order 
to give the reader an idea of some of the possibilities 
suggested in the literature. 

One suggested starting point is the distribution of ef- 
fective mass Mcff = E;^'''''+ET{jl) + ET{j2) + ET{jS) + 
i?T(j4) in the inclusive SUSY sample, which sets the 
approximate mass scale AfsusY = min(mg,TOqi) for the 
strongly interacting sparticles are being produced [[96]. 
and provides a measure of Msusy to 10-15%. 

More detailed information on sparticle masses may 
be accessed by studying specific event^ topologies. For 
instance, the mass of dileptons from Z2 t^l~ Z\ de- 
cays is bounded by m~ — m~ (this bound is even more 

restrictive if Z2 decays via an on-shell slepton)[[H7]. We 
therefore expect an OSSF invariant mass distribution 
to exhibit an edge at m~ — m~ (or below) in any 
sample of SUSY events^ so long as the "spoiler" decay 
modes Z2 — > Z\Z or Z\h are closed. Contamination 
from chargino production can be statistically removed 
by subtracting out the distribution of OSDF dileptons. 
In MHDM models, there may be more than one visi- 
ble mass edge because the Zj, may also be accessible in 
cascade decays. 

In the happy circumstance where production of 
gluinos or a single type of squark is dominant, fol- 
lowed by a string of two-body decays, then further in- 
variant mass edges are possible. One example comes 
from g — > hh\ — > hhZ2 — > hhUZ\\ then one can try to 
combine a 6-jet with the dilepton pair to reconstruct 
the squark-neutralino mass edge: m(b££) < m-j^^ — m~ . 
Next, combining with another 6-jet can yield a gluino- 
neutralino edge: m{bbi£) < rrig — m~ . The reconstruc- 
tion of such a decay chain may be possible as shown 
in Ref. [EB], where other sequences of two-body decays 
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are also examined. In practice, such fortuitous circum- 
stances may not exist, and there are many combinatoric 
issues to overcome as weU. A different study[|98] shows 
that end-point measurements at the LHC will make it 
possible to access the mass difference between the LSP 
and the stau in a mSUGRA scenario where the stau 
co-annihilation mechanism is operative. 

These end-point measurements generally give mass 
differences, not masses. However, by an analysis of 
the decay chain ql — > qZ2 q£^£^ qt^PZi, it 
has been argued[f99] that reconstruction of masses may 
be possible under fortuituous circumstances. More re- 
cently, it has been suggested that it may be possible 
to directly access the gluino and/or squark masses (not 
mass differences) via the introduction of the so-called 
771x2 variable. We will refer the reader to the literature 
for details [[TOD]- 

Mass measurements allow us to check consistency of 
specific SUSY models with a handful of parameters, 
and together with other measurements can readily ex- 
clude such models. But these are not the only interest- 
ing measurements at the LHC. It has been shown that 
if the NLSP of GMSB models decays into a superlight 
gravitino, it may be possible to determine its lifetime, 
and hence the gravitino mass at the LHC[ fTOT] . This 
will then allow one to infer the underlying SUSY break- 
ing scale, a scale at least as important as the weak scale! 
A recent study [[THl] suggests that this is possible even 
when the the decay length of the NLSP is too short 
to be measured. While linear collider experiments will 
ultimately allow the precision measurements that will 
directly determine the new physics to be softly broken 
supersvmmetrv [ 1103] , it will be exciting to analyze real 
LHC data that will soon be available to unravel many of 
the specific details about how (or if) SUSY is actually 
implemented in nature. 

6.8. Measuring DM properties at LHC and ILC 

SUSY discovery will undoubtedly be followed by a 
program (as outlined in Sec. I6.7|l to reconstruct spar- 
ticle properties. What will we be able to say about 
dark matter in light of these measurements? Such 
a study was made by Baltz et al.[ 104J where four 
mSUGRA case study points (one each in the bulk re- 
gion, the HB/FP region, the stau coanihilation region 
and the A-funnel region) were examined for the preci- 
sion with which measurements of sparticle properties 
that could be made at LHC, and also at a -^s = 0.5 
and 1 TeV e+e^ collider. They then adopted a 24- 
parameter version of the MSSM and fit its parame- 
ters to these projected measurements. The model was 
then used to predict several quantities relevant to astro- 



physics and cosmology: the dark matter relic density 

rj— h? , the spin-independent neutralino-nucleon scat- 
^1 ^ 

tering cross section a si (Zip), and the neutralino anni- 
hilation cross section times relative velocity, in the limit 
that V ^ 0: (f7u)|i,_o- The last quantity is the cru- 
cial particle physics input for estimating signal strength 
from neutralino annihilation to anti-matter or gammas 
in the galactic halo. What this yields then is a collider 
measurement of these key dark matter quantities. 

As an illustration, we show in Fig. [11] (taken from 
Ref. [ 1104] ) the precision with which the neutralino 
relic density is constrained by collider measurements for 
the LCC2 point which is in the HB/FP region of the 
mSUGRA model. Measurements at the LHC cannot 
fix the LSP composition, and so unable to resolve the 
degeneracy between a wino-LSP solution (which gives 
a tiny relic density) and the true solution with MHDM. 
Determinations of chargino production cross sections at 
the ILC can easily resolve the difference. It is nonethe- 
less striking that up to this degeneracy ambiguity, ex- 
periments at the LHC can pin down the relic density 
to within ~ 50% (a remarkable result, given that there 
are sensible models where the predicted relic density 
may differ by orders of magnitude!). This improves to 
10-20% if we can combine LHC and ILC measurements. 
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Figure 11. Determination of neutralino relic abundance via 
measurements at the LHC and ILC, taken from Ref. [ \104f - 
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This collider determination of the relic density is very 
important. If it agrees with the cosmological measure- 
ment it would establish that the DM is dominantly ther- 
mal neutralinos from the Big Bang. If the neutralino 
relic density from colliders falls significantly below ([T|), 
it would provide direct evidence for multi-component 
DM- perhaps neutralinos plus axions or other exotica. 
Alternatively, if the collider determination gives a much 
larger value of fi— h^, it could point to a long-lived but 
unstable neutralino and/or non-thermal DM. 

The collider determination of model parameters 
would also pin down the neutralino-nucleon scattering 
cross section. Then if a WIMP signal is actually ob- 
served in DD experiments, one might be able to deter- 
mine the local DM density of neutralinos and aspects of 
their velocity distribution based on the DD signal rate. 
This density should agree with that obtained from as- 
trophysics if the DM in our Galaxy is comprised only 
of neutralinos. 

Finally, a collider determination of {crv)\y^Q would 
eliminate uncertainty on the particle physics side of 
projections for any ID signal from annihilation of neu- 
tralinos in the galactic halo. Thus, the observation of a 
gamma ray and/or anti- matter signal from neutralino 
halo annihilations would facilitate the determination of 
the galactic halo dark matter density distribution. 

7. Some non-SUSY WIMPs at the LHC 

7.1. B}, state from universal extra dimensions 

Models with Universal Extra Dimensions, or UED, 
are interesting constructs which provide a foil for SUSY 
search analyses [[18]. In the 5-D UED theory, one posits 
that the fields of the SM actually live in a 5-D brane 
world. The extra dimension is "universal" since all the 
SM particles propagate in the 5-D bulk. The single 
extra dimension is assumed to be compactified on a 
S1/Z2 orbifold (line segment). After compactification, 
the 4-D effective theory includes the usual SM parti- 
cles, together with an infinite tower of Kaluza-Klein 
(KK) excitations. The masses of the excitations de- 
pend on the radius of the compactified dimension, and 
the first (rt = 1) KK excitations can be taken to be 
of order the weak scale. In these theories, KK-parity 
(— 1)" can be a conserved quantum number. If this 
so-called KK-parity is exact, then the lightest odd KK 
parity state will be stable and can be a DM candidate. 
At tree-level, all the KK excitations in a given level are 
essentially degenerate. Radiative corrections break the 
degeneracy, leaving colored excitations as the heaviest 
excited states and the n = 1 KK excitation of the SM 
U{1)y gauge boson Bj^ as the lightest[ \TU5\ KK odd 



state: in the UED case, therefore, the DM particle has 
spin-1. The splitting caused by the radiative correc- 
tions is also essential to assess how the KK excitations 
decay, and hence are crucial for collider phenomenology 
[[106]. 

The relic density of Bj^ particles has been computed, 
and found to be compatible with observation for certain 
mass ranges of i?^, [ I107j . Also, in UED, the colored ex- 
citations can be produced with large cross sections at 
the LHC, and decay via a cascade to the Bj^ final state. 
Thus, the collider signatures are somewhat reminiscent 
of SUSY, and it is interesting to ask whether it is pos- 
sible to distinguish a jets + leptons + i?™""' signal in 
UED from that in SUSY. Several studies [ [108] answer 
affirmatively, and in fact provide strong motivation for 
the measurement of the spins of the produced new par- 
ticles [ [TUS]. UED DM generally leads to a large rate 
in IceCube, and may also give an observable signal in 
anti-protons and possibly also in photons and positrons [ 
[TS] IllOj . DD is also possible but the SI cross section is 
typically smaller than 10^^ pb. 

7.2. Little Higgs models 

Little Higgs models [ [19] [22] provide an alternative 
method compared to SUSY to evade the quadratic sen- 
sitivity of the scalar Higgs sector to ultra-violet (UV) 
physics. In this framework, the Higgs boson is a pseudo- 
Goldstone boson of a spontaneously broken global sym- 
metry that is not completely broken by any one cou- 
pling, but is broken when all couplings are included. 
This then implies that there is quadratic sensitivity to 
UV physics, but only at the multi-loop level. Specific 
models where the quadratic sensitivity enters at the 
two-loop level should, therefore, be regarded as low en- 
ergy effective theories valid up to a scale A ~ 10 TeV, 
at which a currently unknown, and perhaps strongly- 
coupled UV completion of the theory is assumed to 
exist. Models that realize this idea require new TeV- 
scale degrees of freedom that can be searched for at 
the LHC: new gauge bosons, a heavy top-like quark, 
and new spin-zero particles, all with couplings to the 
SM. These models, however, run into phenomenologi- 
cal difficulties with precision EW constraints, unless a 
discrete symmetry- dubbed r-parity[[20j- is included. 
SM particles are then T-even, while the new particles 
are T-odd. 

We will set aside the issue (mentioned earlier) of 
whether T-parity conservation is violated by anoma- 
lies [[23], and assume that a conserved T-parity can be 
introduced['2T. In this case, the lightest T-odd particle 
Ah - the Little Higgs partner of the hypercharge gauge 
boson with a small admixture of the neutral W^h boson 
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- is stable and yields the observed amount of DM for 
a reasonable range of model parameters[ IllOj . In this 
case, the DM particle has spin-1, though other cases 
with either a spin-^ or spin-0 heavy particle may also 
be possible. Ah can either annihilate with itself into 
vector boson pairs or tt pairs via s-channcl Higgs ex- 
change, or into top pairs via exchange of the heavy T- 
odd quark in the t-channel. Co-annihilation may also 
be possible if the heavy quark and Ah are sufficiently 
close in mass. Signals at the LHC[ llllj mainly come 
from pair production of heavy quarks, and from single 
production of the heavy quark in association with Ah- 
These lead to low jet multiplicity events plus E^^^^. 
The E™^^^ comes from the escaping Ah particle, which 
must be the endpoint of all T-odd particle decays If 
Ah is the dominant component of galactic DM, we will 
generally expect small DD and ID rates for much the 
same reasons that the signals from the bino LSP tend 
to be small [ [TTO] : see, however, Ref.[ lll2] for a different 
model with large direct detection rate. 

8. Outlook 

The union of particle physics, astrophysics and cos- 
mology has reached an unprecedented stage. Today we 
are certain that the bulk of the matter in the universe 
is non- luminous, not made of any of the known par- 
ticles, but instead made of one or more new physics 
particles that do not appear in the SM. And though we 
know just how much of this unknown dark matter there 
is, we have no idea what it is. Today, many theoreti- 
cal speculations which seek to answer one of the most 
pressing particle physics puzzles, "What is the origin of 
EWSB and how can we embed this into a unified the- 
ory of particle interactions?" automatically also point 
to a resolution of this 75 year old puzzle as to what the 
dominant matter component of our universe might be. 
Particle physicists have made many provocative sug- 
gestions for the origin of DM, including supersymmetry 
and extra spatial dimensions, ideas that will completely 
change the scientific paradigm if they prove to be right. 

The exciting thing is that many of these speculations 
will be directly tested by a variety of particle physics 
experiments along with astrophysical and cosmologi- 
cal searches. The Large Hadron Collider, scheduled to 
commence operation in 2008, will directly study parti- 
cle interactions at a scale of 1 TeV where new matter 



''We note here that it is also possible to construct so-called twin- 
Higgs modelsf 11131 where the Higgs sector is stabilized via new 
particles that couple to the SM Higgs doublet, but are singlets 
under the SM gauge group. In this case, there would be no 
obvious new physics signals at the LHC. 



states are anticated to exist for sound theoretical rea- 
sons. These new states may well be connected the DM 
sector, and so in this way the LHC can make crucial 
contributions to not only particle physics, but also to 
cosmology. 

Any discovery at LHC of new particles at the TeV 
scale will make a compelling case for the construction 
of a lepton collider to study the properties of these par- 
ticles in detail and to elucidate the underlying physics. 
Complementary to the LHC, there are a variety of 
searches for signals from relic dark matter particles ei- 
ther locally or dispersed throughout the galactic halo. 
The truly unprecedented thing about this program is 
that if our ideas connecting DM and the question of 
EWSB are correct, measurements of the properties of 
new particles produced at the LHC (possibly comple- 
mented by measurements at an electron-positron linear 
collider) may allow us to independently infer just how 
much DM there is in the universe, and quantitatively 
predict what other searches for DM should find0 

Particle physics, cosmology and astrophysics are 
rapidly obliterating their boundaries and merging into 
a single discipline. The ACDM model that has emerged 
posits that 70% of the energy budget of the Universe 
is contained in so-called dark energy, weird stuff with 
negative pressure that is completely different from any- 
thing that we have ever encountered! Thus, not only 
are the particles we are made of a small fraction of the 
total matter content of the Universe, most of the en- 
ergy of the universe appears to be in non-material dark- 
energy, extending even further the Copernican princi- 
pleO This ACDM framework is being incisively probed 
by observation, and may possibly need modification. 
The nature of dark energy is a completely open ques- 
tion. Experiments over the next decade or two will, 
we expect, reveal the identity of dark matter and, we 
hope, will provide clues as to the origin of dark en- 
ergy. This unprecendented synthesis of the physics of 
both the largest and smallest scales observable in na- 
ture should make the next twenty years very exciting! 
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■^These studies have only just begun, and have only been car- 
ried out in the context of supersymetry, which unlike extra- 
dimensional or Little Higgs models, is a complete theory, valid 
up to very high energy. 

^Our colleagues who subscribe to the multiverse view carry this 
yet further, suggesting that our Universe is just one of many. 
Unlike for the ideas discussed here, we are not aware of possible 
tests for this view. 



Dark matter and the LHC 



23 



REFERENCES 

1. For reviews, see e.g. C. Jungman, M. Kamionkowski 
and K. Gnest,Phy.s. Kept. 267 (1996) 195; A. La- 
hanas, N. Mavromatos and D. Nanopoulos, Int. 
J. Mod. Phys. D 12 (2003) 1529; M. Drees, hep- 
ph/0410113; K. Olive, "Tasi Lectures on Astropar- 
ticle Physics", astro-ph/0503065; G. Bertone, D. 
Hooper and J. Silk, Phys. Kept. 405 (2005) 279. 

2. F. Zwicky, Helvetica Physica Acta 6 (1933) 110; see 
also Astrophys. J. 86 (1937) 217. 

3. V. Rubin and W. K. Ford, Astrophys. J. 159 (1970) 
359; V. Rubin, N. Thonnard and W. K. Ford, As- 
trophys. J. 238 (1980) 471. 

4. A. Conley et al. Astrophys. J. 644 (2006) 1. 

5. D. N. Spergel et al. (WMAP Collaboration), Astro- 
phys. J. Supp., 170 (2007) 377. 

6. B. Lee and S. Weinberg, Phys. Rev. Lett. 39 (1977) 
165. 

7. Report of the DMSAG panel on 
http : / /www. science . doe . gov /hep /hepapjreports. shtm 
See also, L. Roszkowski, Pramama, 62 (2004) 389. 

8. A. Kusenko, Phys. Rev. Lett. 97 (2006) 241301. 

9. P. Chen, Mod. Phys. Lett.A 19 (2004) 1047. 

10. S. Coleman, Nucl. Phys. B 262 (1985) 263, 
[Erratum-ibid. B269 (1986) 744. 

11. A. Kusenko and M. Shaposhnikov, Phys. Lett. B 
418 (1998) 46; A. Kusenko, V. Kuzmin and M. 
Shaposhnikov, Phys. Rev. Lett. 80 (1998) 3185. 

12. E. Kolb, D. Chung and A. Riotto, Phys. Rev. D 59 
(1999) 023501 and hcp-ph/9810361. 

13. R. Pcccei and H. Quinn, Phys. Rev. Lett. 38 (1977) 
1440 and Phys. Rev. D 16 (1977) 1791; S. Wein- 
berg, Phys. Rev. Lett. 40 (1978) 223; F. Wilczek, 
Phys. Rev. Lett. 40 (1978) 279. 

14. For a recent review, see e.g. S. J. Asztalos, L. Rosen- 
berg, K. van Bibber, P. Sikivie and K. Zioutas, Ann. 
Rev. Nucl. Part. Sci.bQ (2006) 293. 

15. For recent results, see L. Duffy et al., Phys. Rev. 
Lett. 95 (2005) 091304. 

16. H. Goldberg, Phys. Rev. Lett. 50 (1983) 1419; J. El- 
lis, J. Hagclin, D. Nanopoulos and M. Srednicki, 
Phys. Lett. B 127 (1983) 233; J. Ellis, J. Hagelin, 
D. Nanopoulos, K. Olive and M. Srednicki, Nucl. 
Phys. B 238 (1984) 453. 

17. G. Servant and T. Tait, New Jou. Phys.4 (2002) 99. 

18. For a review, see D. Hooper and S. Profumo, Phys. 
Rept. 453 (2007) 29. 

19. N. Arkani-Hamed, A. Cohen and H. Geoigi, Phys. 
Lett. B 513 (2001) 232; N. Arkani-Hamed, A. Co- 
hen, E. Katz and A. Nelson, J. High Energy Phys. 
07 (2002) 034. 



20. H. C. Cheng and L Low, J. High Energy Phys. 0309 
(2003) 051. 

21. J. Hubisz and P. Meade, Phys. Rev. D 71 (2005) 
035016. 

22. For reviews, see M. Schmaltz, Ann. Rev. Nucl. Part. 
Sci.55 (2005) 229 and M. Perelstein, Prog. Part. 
Nucl. Phys.58 (2007) 247. 

23. C. HiU and R. HiU Phys. Rev. D 75 (2007) 115009. 

24. H.-C. Cheng, arXiv:07 10.3407 [hep-ph]. 

25. J. Feng and J. Kumar. [arXiv:0803.4196l [hep-ph]. 

26. K. Jcdamzik, M. LeMoine and G. Moultaka, 
JCAP0607 (2006) 010. 

27. S. Heinemeyer, W. Holik, D. Stockinger, A. Weber 
and G. Weiglein, J. High Energy Phys. 0608 (2006) 
052. 

28. H. Baer and X. Tata, Weak Scale Super symmetry: 
From Superfields to Scattering Events, (Cambridge 
University Press, 2006). 

29. M. Drees, R. Godbole and P. Roy, Sparticles, 
(World Scientific, 2004). 

30. P. Binetruy, Supersymmetry (Oxford, 2006). 

31. S. Dimopoulos and H. Georgi, Nucl. Phys. B 193 
(1981) 150. 

32. B.C. AUanach, S. Kraml and W. Porod, J. High 
Energy Phys. 03 (2003) 016; G. Belanger, S. Kraml 
and A. Pukhov, Phys. Rev. D 72 (2005) 015003; S. 
Kraml and S. Sekmen in: Physics at TeV Collid- 
ers 2007, BSM working group report, in prep.; see 
http : / / cern . ch/kraml/ comparison. 

33. M. Dine, A. Nelson, Y. Nir and Y. Shirman, Phys. 
Rev. D 53 (1996) 2658; for a review, see G. Giudice 
and R. Rattazzi, Phys. Rept. 322 (1999) 419. 

34. L. Randall and R. Sundrum, Nucl. Phys. B 557 
(1999) 79; G. Giudice, M. Luty, H. Murayama and 
R. Rattazzi, J. High Energy Phys. 9812 (1998) 027. 

35. S. Kachru, R. KaUosh, A. Linde and S. P. Trivedi, 
Phys. Rev. D 68 (2003) 046005; K. Choi, A. 
Falkowski, H. P. Nilles, M. Olechowski and S. Poko- 
rski, J. High Energy Phys. 0411 (2004) 076; K. 
Choi, A. Falkowski, H. P. Nilles and M. Olechowski, 
Nucl. Phys. B 718 (2005) 113; K. Choi, K-S. 
Jeong and K. Okumura, J. High Energy Phys. 0509 
(2005) 039; H. Baer, E. Park, X. Tata and T. Wang, 
J. High Energy Phys. 0706 (2007) 033, and refer- 
ences therein. 

36. L. Everett, I.-W. Kim and K. Zurek, 
arXiv:0804:0592. 

37. P. Gondolo, J. Edsjo, P. UUio, L. Bergstrom, 
M. Schelke and E. A. Baltz, JCAP 0407 (2004) 
008. 

38. G. Belanger, F. Boudjema, A. Pukhov and A. Se- 
menov, Comput. Phys. Commun.\7A (2006) 577; 



24 



H. Baer and X. Tata 



Comput. Phys. Commun.176 (2007) 367. 

39. IsaRED, by H. Baer, C. Balazs and A. Belyaev, J. 
High Energy Phys. 0203 (2002) 042. 

40. ISAJET, by H. Baer, F. Paige, S. Protopopescu and 
X. Tata, hep-ph/0312045; see also H. Baer, J. Fer- 
randis, S. Kraml and W. Porod, Phys. Rev. D 73 
(2006) 015010. 

41. H. Baer and M. Brhlik, Phys. Rev. D 53 (1996) 
597; V. Barger and C. Kao, Phys. Rev. D 57 (1998) 
3131. 

42. J. Ellis, T. Falk and K. Olive, Phys. Lett. B 444 
(1998) 367; J. Ellis, T. Falk, K. Olive and M. Sred- 
nicki, Astropart. Phys. 13 (2000) 181; M.E. Gomez, 
G. Lazarides and C. Pallis, Phys. Rev. D 61 (2000) 
123512 and Phys. Lett. B 487 (2000) 313; A. La- 
hanas, D. V. Nanopoulos and V. Spanos, Phys. Rev. 
D 62 (2000) 023515; R. Arnowitt, B. Dutta and 
Y. Santoso, Nucl. Phys. B 606 (2001) 59; see also 
Ref. [|39]. 

43. K. L. Chan, U. Chattopadhyay and P. Nath, Phys. 
Rev. D 58 (1998) 096004; J. Feng, K. Matchev 
and T. Moroi, Phys. Rev. Lett. 84 (2000) 2322 and 
Phys. Rev. D 61 (2000) 075005; see also H. Baer, 
C. H. Chen, F. Paige and X. Tata, Phys. Rev. 
D 52 (1995) 2746 and Phys. Rev. D 53 (1996) 
6241; H. Baer, C. H. Chen, M. Drees, F. Paige 
and X. Tata, Phys. Rev. D 59 (1999) 055014; 
for a model-independent approach, see H. Baer, 
T. Krupovnickas, S. Profumo and P. UUio, J. High 
Energy Phys. 0510 (2005) 020. 

44. M. Drees and M. Nojiri, Phys. Rev. D 47 (1993) 
376; H. Baer and M. Brhlik, Phys. Rev. D 57 
(1998) 567; H. Baer, M. Brhlik, M. Diaz, J. Fer- 
randis, P. Mercadante, P. Quintana and X. Tata, 
Phys. Rev. D 63 (2001) 015007; J. Ellis, T. Falk, 
G. Ganis, K. Olive and M. Srednicki, Phys. Lett. B 
510 (2001) 236; L. Roszkowski, R. Ruiz de Austri 
and T. Nihei, J. High Energy Phys. 0108 (2001) 
024; A. Djouadi, M. Drees and J. L. Kneur, J. 
High Energy Phys. 0108 (2001) 055; A. Lahanas 
and V. Spanos, Eur. Phys. J. C 23 (2002) 185. 

45. R. Arnowitt and P. Nath, Phys. Rev. Lett. 70 
(1993) 3696; H. Baer and M. Brhlik, Ref. [ |4T]; 
A. Djouadi, M. Drees and J. Kneur, Phys. Lett. 
B 624 (2005) 60. 

46. C. Bohni, A. Djouadi and M. Drees, Phys. Rev. 
D 30 (2000) 035012; J. R. Ellis, K. A. Olive 
and Y. Santoso, Astropart. Phys. 18 (2003) 395; 
J. Edsjo, et at, JCAP 0304 (2003) 001. 

47. H. Baer, A. Mustafayev, E. Park and X. Tata, 
JCAP0701, 017 (2007). 

48. H. Baer, A. Mustafayev, E. Park and X. Tata, 



arXiv:0802:3384. 

49. Z. Ahmed et o/., arXiv:0802:3530 [astro-ph]. 

50. J. Angle et al, Phys. Rev. Lett. 100 (2008) 021303. 

51. E. Behnk et al, arXiv:0804:2886 [astro-ph] 

52. D. Feldman, Z. Liu and P. Nath, Phys. Lett. B 662 
(2008) 190. 

53. H. Baer, C. Balazs, A. Belyaev and J. OTarrill, 
JCAP 0309, (2003) 007. 

54. R. Schnee, (CDMS Collaboration); A. M. Green, 
JCAP 0708 (2007) 022; C-L. Shan and M. Drees, 
arXiv:0710.4296 [hep-ph]. 

55. H. Baer, A. Belyaev, T. Krupovnickas and 
J. O'FarriU, JCAP 0408 (2004) 005. 

56. H. Baer and J. O'FarriU, JCAP0404, 005 (2004). 

57. P. Sreekumar et al. [EGR ET Collaboration ] , Astro - 
phys. J. 494, 523 (1998) |arXiv:astro-ph/970925';i . 

58. W. de Boer, M. Herold, C. Sander, V. Zhukov, 
A. V. Gladyshev and D. I. Kazakov, 
|arXiv:astro-ph/ 0408272, 

59. F. W. Stecker, S. D. Hunter and D. A. Kniffcn, 
arXiv:0705.4311 [astro-ph]. 

60. H. Baer, A. Belyaev and H. Sumniy, 
arrXiVr0802:3r27 [hep-ph]. 

61. S. Weinberg, Phys. Rev. Lett. 48 (1982) 1303; R. H. 
Cyburt, J. Ellis, B. D. Fields and K. A. Olive, Phys. 
Rev. D 67 (2003) 103521; K. Jedamzik, Phys. Rev. 
D 70 (2004) 063524; M. Kawasaki, K. Kohri and 
T. Moroi, Phys. Lett. B 625 (2005) 7 and Phys. 
Rev. D 71 (2005) 083502. K. Kohri, T. Moroi and 
A. Yotsuyanagi, Phys. Rev. D 73 (2006) 123511. 

62. H. Pagels and J. Primack, Phys. Rev. Lett. 
48 (1982) 223; J. Feng, A. Rajaraman and 
F. Takayama, Phys. Rev. Lett. 91 (2003) 011302 
and Phys. Rev. D 68 (2003) 085018. 

63. M. Bolz, A. Brandenburg and W. BuchmuUer, 
Nucl. Phys. B 606 (2001) 518; J. Pradler and F. 
Steffen, hep-ph/0608344. 

64. J. Feng, S. Su and F. Takayama, Phys. Rev. D 70 
(2004) 075019. 

65. W. BuchmuUer, P. Di Bari and M. Plumacher, An- 
na! Phys. 315 (2005) 305. 

66. W. BuchmuUer, L. Covi, J. Kersten, K. Schmidt- 
Hoberg, JCAP0611 (2006) 007; W. BuchmuUer, L. 
Covi, K. Hamaguchi, A. Ibarra and T. Yanagida, 
J. High Energy Phys. 0703 (2007) 037. 

67. J. E. Kim and H. P. Miles, 

68. L. Covi, J. E. Kim and L. Roszkowski, Phys. Rev. 
Lett. 82 (1999) 4180; L. Covi, H. B. Kim, J. E. 
Kim and L. Roszkowski, J. High Energy Phys. 0105 
(2001) 033. 

69. A. Brandenburg and F. Steffen, JCAP0408 (2004) 
008. 



Dark matter and the LHC 



25 



70. G. Lazarides and Q. Shafi, Phys. Lett. B 258 (1991) 
305; K. Kumekawa, T. Moroi and T. Yanagida, 
Prog. Theor. Phys. 92 (1994) 437; T. Asaka, K. 
Hamaguchi, M. Kawasaki and T. Yanagida, Phys. 
Lett. B 464 (1999) 12. 

71. H. Baer and H. Summv IarXlv:0803.0510 l (2008). 

72. H. Baer, S. Kraml, S. Sekmen and H. Summy, 
JHEP0803, 056 (2008). 

73. H. Baer and X. Tata, Phys. Lett. B 160 (1985) 159. 

74. H. Baer, J. Ellis, G. Gelmini, D. V. Nanopoulos 
and X. Tata, Phys. Lett. B 161 (1985) 175; G. 
Gamberini, Z. Physik C 30 (1986) 605; H. Baer, 
V. Barger, D. Karatas and X. Tata, Phys. Rev. D 
36 (1987) 96; H. Baer, X. Tata and J. Woodside, 
Phys. Rev. D 45 (1992) 142. 

75. H. Baer, D. Dzialo-Karatas and X. Tata, Phys. Rev. 
D 42 (1990) 2259. 

76. H. Baer, C. H. Chen, F. Paige and X. Tata, Phys. 
Rev. D 50 (1994) 4508. 

77. H. Baer, C. H. Chen, F. Paige and X. Tata, Phys. 
Rev. D 49 (1994) 3283 

78. H. E. Haber and D. Wyler, Nucl. Phys. B 323 
(1989) 267; S. Ambrosanio and B. Mele, Phys. Rev. 
D 53 (1996) 2541 and Phys. Rev. D 55 (1997) 1399 
[Erratum-ibid. D56, 3157 (1997)]; H. Baer and T. 
Krupovnickas, J. High Energy Phys. 0209 (2002) 
038. 

79. H. Baer and X. Tata, Phys. Rev. D 47 (1993) 2739. 

80. H. Baer, C. Chen, M. Drees, F. Paige and X. Tata, 
Phys. Rev. Lett. 79 (1997) 986. 

81. H. Baer, C. Chen, M. Drees, F. Paige and X. Tata, 
Phys. Rev. D 59 (1999) 015010. 

82. H. Baer, M. Bisset, X. Tata and J. Woodside, Phys. 
Rev. D 46 (1992) 303. 

83. H. Baer, M. Bisset, D. Dicus, C. Kao and X. Tata, 
Phys. Rev. D 47 (1993) 1062; H. Baer, M. Bisset, 
C. Kao and X. Tata, Phys. Rev. D 50 (1994) 316. 

84. M. Muhlleitner, A. Djouadi and Y. Mambrini, 
Comput. Phys. Commun.WS (2005) 46. 

85. W. Porod, Comput. Phys. Commun.153 (2003) 
275. 

86. T. Sjostrand, S. Mrenna and P. Skands, J. High 
Energy Phys. 0605 (2006) 026. 

87. G. Corcella et ai, J. High Energy Phys. 0101 
(2001) 010. 

88. F. Maltoni and T. Stelzer, J. High Energy Phys. 
0302 (2003) 027; J. AlwaU et ai, J. High Energy 
Phys. 0709 (2007) 028. 

89. A. Pukhov et al, hep-ph/9908288. 

90. W. Kilian, T. Ohl and J.Reuter, 'arXiv:0708.4233l 

91. J. AlwaU et ai, Comput. Phys. Commun.176 
(2007) 300. 



92. H. Baer, X. Tata and J. Woodside, Ref.[7^. 

93. P. Mercadante, J. K. Mizukoshi and X. Tata, 
Phys. Rev. D 72 (2005) 035009; S. P. Das 
et al. arXiv:0708:2048 [hep-ph]; R. Kadala, 
P. Mercadante, J. K. Mizukoshi and X. Tata, 
arXiv:0803:0001 [hep-ph]. 

94. H. Baer, A. Belyaev, T. Krupovnickas and X. Tata, 
J. High Energy Phys. 0402 (2004) 007; H. Baer, 
T. Krupovnickas and X. Tata, J. High Energy Phys. 
0406 (2004) 061. 

95. H. Baer, H. Prosper and H. Summy, Phys. Rev. D 
77 (2008) 055017. 

96. L Hinchhffe et al, Phys. Rev. D 55 (1997) 5520 
and Phys. Rev. D 60 (1999) 095002. 

97. H. Baer, K. Hagiwara and X. Tata, Phys. Rev. D 
35 (1987) 1598; H. Baer, D. Dzialo-Karatas and X. 
Tata, Phys. Rev. D 42 (1990) 2259; H. Baer, C. 
Kao and X. Tata, Phys. Rev. D 48 (1993) 5175; H. 
Baer, C. H. Chen, F. Paige and X. Tata, Phys. Rev. 
D 50 (1994) 4508. 

98. R. Arnowitt et al. Phys. Lett. B 639 (2006) 46 and 
Phys. Lett. B 649 (2007) 73. 

99. H. Bachacou, I. HinchlifFe and F. Paige, Phys. Rev. 
D 62 (2000) 015009; Atlas CoUaboration, LHCC 
99-14/15. 

100. C. Lester and D Summers, Phys. Lett. B 463 (1999) 
99; A. Barr, C. Lester and P. Stephens, J. Phys. bf 
G29 (2003) 2343; C. Lester and A. Barr, J. High 
Energy Phys. 0712 (2007) 102; W. Cho, K. Choi, 
Y. Kim and C. Park, arXiv:0711.4526, 

101 K. Kawagoe et al. Phys. Rev. D 69 (2004) 035003; 
S. Ambrosanio et al. J. High Energy Phys. 0101 
(2001) 014. 

102H. Hamaguchi, S. Shirai and T. Yanagida, 

IarXiv:0712.2462, 
103J. Feng et al. Phys. Rev. D 52 (1995) 1418; M. No- 

jiri, K. Fujii and T. Tsukamoto, Phys. Rev. D 54 

(1996) 6756. 

104E. Bahz, M. Battagha, M. Peskin and T. Wiz- 
ansky, Phys. Rev. D 74 (2006) 103521. See also 
R. Arnowitt et al, arXiv: 0802. 2968 [hep-ph] for a 
similar study in the stau co-annihilation region. 

105B. C. Cheng, K. Matchev and M. Schmahz, Phys. 
Rev. D 66 (2002) 036005. 

106.T. Rizzo, Phys. Rev. D 64 (2001) 095010; C. 
Macescanu, C. McMuUen and S. Nandi, Phys. Rev. 
D 66 (2002) 015009. 

107.See Ref. [\l7\; K. Kong and K. Matchev, J. High 
Energy Phys. 0601 (2006) 038. 

108H. C. Cheng, K. Matchev and M. Schmahz, Phys. 
Rev. D 66 (2002) 056006; A. Datta, K. Kong and 
K. Matchev, Phys. Rev. D 72 (2005) 096006. 



26 



H. Baer and X. Tata 



109A. Alvcs, O. Eboli and T. Plehn, Phys. Rev. D 74 
(2006) 095010. 

IIOA. Birkedal, A Noble, M. Perelstein and A. Spray, 
Phys. Rev. D 74 (2006) 035002; D. Hooper and 
G. Zaharijas, Phys. Rev. D 75 (2007) 035010. 

111. C.-S. Chen, K. Cheung and T. C. Yuan, Phys. 
Lett. B 664 (2007) 158; T. Han, H. Logan and L.- 
T. Wang, J. High Energy Phys. 0601 (2006) 099. 

112. Y. Bai, arXiv:0801:1662. 

113. Z. Chacko, H-S. Goh and R. Harnik, Phys. Rev. 
Lett. 96 (2006) 231802; see E. DoUe and S. Su, 
Phys. Rev. D 77 (2008) 075013 for an analysis of 
DM in a twin-Higgs scenario. 



